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Abstract 

A  steady  state  heat  transfer  technique  is  developed  and  evaluated  for  detecting 
boundary  layer  transition  on  a  flat  plate  in  incompressible  flow.  The  method  involves 
adhering  eiK^Mulaied  temperature  sensitive  liquid  crystals  to  a  constant  heat  flux 
surface.  A  huter  conqrosed  of  unidirectioiud  carbon  fibres  is  devdoped  and  tested 
with  dte  aim  of  in-flight  boundary  layer  transition  detection  <»)  a  natural  laminar  flow 
nacelle.  The  individual  and  combirted  influences  of  surface  heating  and  favourable 
pressure  gradient  cm  boundary  layer  transititni  are  considered.  Heating  is  found  to  be 
a  destabilising  influence  on  die  boundary  layer  while  a  favourable  pressure  gradient 
is  a  stabilising  influence.  A  Stanton  number  correlates  die  movement  of  transition 
with  heating  for  zero  pressure  gradient  flat  plate  boundary  layers.  No  similar 
correlation  is  found  when  a  favourable  pressure  gradient  acccunpanies  die  heating. 
Heating  is  more  destabilising  in  a  favourable  pressure  gradient  than  in  zero  pressure 
gradient.  Shear  sensitive  liquid  crystal  is  used  to  detect  transition  and  to  obtain,  for 
the  first  time,  quantitative  measutemaits  of  surface  shear  stress  in  a  transitional 
boundary  layer.  This  involves  calibrating  the  time  required  for  the  liquid  crystal  film 
to  experience  a  texture  change  when  subjected  to  a  shearing  force.  Using  this 
technique,  shear  sensitive  liquid  crystal  shows  pranise  for  full  coverage  measurement 
of  surface  shear  stress  in  low  speed  flows. 
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Nomenclature 

' 

C, 

specific  heat  at  constant  piessuie,  J^-K 

h 

coefficient  of  heat  transfer, 

i 

electric  current,  Anoqps 

k 

thermal  conductivity,  W/m  K 

K 

acceleration  parameter 

P 

piessttie, 

Pr 

Prandtl  number 

<l 

heat  flux,  W/m^ 

R 

electrical  resistance,  £2 

Re 
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St 
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t 
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t 

time,  sec 

T 
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n 

liquid  crystal  colour  play  temperature,  K 

T, 

gas  temperature,  K 

T. 

surface  teinperatue,  K 

Tu 

turbulence  intensity 

U 

fluid  velocity,  m/sec 

u 

fluid  velocity  conqwnent,  m/sec 

“/ 

local  fluid  velocity  component,  m^ec 

JC 

distance  parallel  to  surface,  m 
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distance  noimal  to  surface,  m 
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temperature  coefficient  of  resistance,  f2/K 

Y 
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boundary  layer  thickness,  m 
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Subscripts 


ad 

adiabatic 

cond 

conduction 

conv 

convection 

crit 

critical 

ref 

reference  conditions 

t 

transition  start 

T 

transition  end 

w 

conditions  at  the  wall 

Chapter  1:  Introduction  and  Literature  Survey 


1.1  Natural  Lananar  Flow 

The  Oxford  University  Engineering  Laboratory  is  to  be  involved  in  experiments 
to  detect  boundary  layer  transition  in  flight  tests  on  a  natural  laminar  flow  (NLF) 
engine  nacelle.  Tte  skin  friction  drag  associated  widi  a  laminar  boundary  layer  is 
substantially  smaller  than  that  of  an  equivalent  length  of  turbulent  boundary  layer. 
Laminar  flow  technology  is,  therefore,  directed  at  reducing  viscous  drag  on  aircraft 
by  delaying  the  point  of  boundary  layer  transition  from  laminar  to  turbulent  flow. 
Transition  may  be  delayed  by  using  either  active  or  passive  means.  Ihe  latter  is 
advantageous  in  that  the  power  requirements  of  active  systems,  such  as  suction, 
diminish  the  benefits  of  delayed  transition.  One  such  passive  measure  is  natural 
laminar  flow  technology,  in  which  die  geometry  of  an  aerodynamic  body  is  designed 
so  as  to  delay  transition. 

It  has  been  established  that  die  point  of  transition  in  the  boundary  layer  is 
strongly  influenced  by  the  pressure  gradient  in  the  external  streariL  A  negative  (or 
favourable)  pressure  gradient  in  die  stream-wise  direction  leads  to  delayed  transititm 
and  an  adverse  pressure  gradient  advances  transititm  location.  NLF  technology  seeks 
to  extend  the  length  of  laminar  flow  by  developing  a  body  geometry  which  has  a 
favourable  pressure  gradient  over  as  much  of  its  length  as  possible  (Fig  1.1).  There 
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is,  however,  a  trade-off  to  be  considered  when  attenq>ting  to  reduce  drag  by  this 
means.  A  favourable  pressure  gradient  is  formed  at  the  expense  of  an  increase  in 
wing  thickness  (or  nacelle  diameter)  which  leads  to  increased  form  drag. 

HASTINGS  et  al.  (1987)  have  noted  that  the  external  nacelle  on  laige 
turbofan-powered  aircraft  are  particularly  attractive  candidates  for  the  application  of 
natural  laminar  flow  technology.  These  nacelles  have  lower  characteiistic  Rejmolds 
numbers  than  other  possible  candidates,  such  as  fuselages  and  wings.  Furthermore, 
nacelles  are  not  required  to  provide  lift  and  can,  therefore,  be  shaped  to  have  pressure 
gradients  favourable  to  laminar  boundary  layers  without  the  concern  for  lift  and 
moment  characteristics  which  influences  the  design  of  NLF  wings. 

The  practical  importance  of  NLF  woik  is  considerable.  By  contouring  the  NLF 
nacelle  so  as  to  have  accelerating  flow  over  as  much  as  70%  of  its  length,  HASTINGS 
et  al.  suggest  that  nacelle  ffiction  drag  may  be  reduced  by  as  much  as  40-30%.  This 
equates  to  a  reduction  in  total  aircraft  drag  and  cruise  fuel  bum  of  1-2%  for  a  large 
commercial  transport.  The  potential  savings  in  fuel  costs  (see  Fig.  1.2)  as  well  as 
benefits  in  terms  of  increased  flight  range  and  payload  and  reduced  environiriental 
pollution  are  substantial.  'These  benefits  may  be  increased  considerably  by  NLF 
triplications  to  other  aircraft  structures  (Fig.  1.3).  GODARD  (1989)  has  noted  that 
viscous  drag  accounts  for  -  30%  of  total  aircraft  drag.  SZODRUCH  (1991)  estinutes 
that  NLF  technology  offers  a  potential  reduction  of  this  drag  by  13.3%. 


Although  there  is  presently  no  commercial  aircraft  desigiwd  for  laminar 
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boundary  layers,  there  have  been  a  number  of  laminar  flow  projects  over  the  last 
decade.  Since  the  early  1980’s,  General  Electric  and  NA.SA’s  Langley  Research 
Center  have  combined  to  peifoim  NLF  wind  nmnel  and  flight  tests.  Tests  were 
performed  at  unit  Reynolds  numbers  of  -  3.5-6.0  x  10*  m  '  on  Citation  and  OV-IB 
aircraft.  Phase  I  of  these  experiments  suggests  that  NLF  could  be  maintained  as  far 
aft  as  S0%  of  nacelle  length. 

COUSTEEX  et  al.  (1991)  have  reported  that  ONERA/CERT  have  performed 
laminar  flow  flight  experiments  on  the  vertical  tail  of  a  Falcon  SO  aircraft.  Flight 
conditions  included  the  transonic  range  with  unit  Reynolds  numbers  of  4-7  x  10*  m  '. 


In  198S  Germany  also  initiated  a  national  lesearch  program  intended  to  develc^ 
transonic  laminar  flow  techirolog;;.  REDEKER  et  al.  (1990),  HORSTMANN  et  al. 
(1990  and  1990a)  and  HENKE  et  al.  (1990)  describe  flight  tests  employing  a  wing 
glove  on  a  light  aircraft  and  the  corresponding  wiird  tunrrel  tests  using  a  1/2  scale 
model  of  the  wing.  The  program  has  thus  far  been  focused  on  performing  flight  and 
wind  tunnel  transition  investigations  on  wings.  The  goal  is  to  use  die  results  to 
develop  a  reliable  transition  prediction  procedure  and  a  leasonaUe  correlatitm  between 
flight  and  wind  tuimel  data.  Tests  were  conducted  at  altitudes  of  4000-10000  ft,  and 
at  Reynolds  numbers  based  <m  mean  aerodynamic  chord  lengdi  of  3-10  x  l(f.  This 
program  has  demonstrated  the  NLF  concept  with  laminar  flow  in  flight  of  up  to  30% 
of  the  chord. 


I 


The  experiments  presented  herein  are  related  to  tests  which  are  to  be  petfonned 
on  an  NLF  nacelle  as  pait  of  a  joint  ventnre  involving  Rolls-Royce  pic.,  DLR,  Hurel- 
Ehibois  (UK)  ltd.,  Oxford  University,  and  MTU.  Two  series  of  flight  tests  are  to  be 
performed  as  a  means  of  demonstrating  NLF  technology.  In  die  first  tests, 
measurements  will  be  made  of  drag  and  transition  locatitm  on  the  nacelle  of  a  Rolls- 
Royce  M45  turbofan  engine  over  a  range  of  flight  conditions.  This  nacelle  is 
constiucted  to  the  original  coordinates  for  the  machine  and  to  typical  standards  for  a 
modem  production  engine.  It  is  referred  to  as  the  datum  nacelle.  This  nacelle  will 
then  be  replaced  with  a  second  nacelle  whose  geometry  has  been  designed  by  Rolls- 
Royce  and  DLR  to  extend  laminar  flow.  Transition  location  and  drag  wiU  be 
measured  on  this  NLF  nacelle  over  a  range  of  flight  conditions  and  engine  settings 
similar  to  those  for  which  the  datum  nacelle  was  tested.  This  program  will  enable 
direa  con:q)arisons  to  be  made  between  the  datum  and  NLF  nacelles  under  similar 
free-stream  conditions  and  power  settings  and  with  identical  engines.  DLR’s  VFW614 
(Fig.  1.4)  will  serve  as  the  test  vehicle  and  the  port-side  M45  engine  will  be  instru¬ 
mented. 

DLR  and  Oxford  University  are  charged  with  transition  detection.  DLR  will 
utilise  infiaied  imaging  while  Oxford  will  measure  transition  using  temperature 
sensitive  liquid  ciystals  tqtplied  to  a  carbon-fibre  heater  vriiich  is  integrated  into  the 
structure  of  the  nacelle.  Resistance  tenqterature  detectors  beneath  the  carbon-fibre 
beater  will  provide  an  alternative  means  of  measuring  surface  temperature  and  serve 
as  feedback  units  to  the  heater’s  temperature  control  system.  A  series  of  hot  film 
gauges  develtqted  and  manufactured  by  Oxford  University  will  provide  a  third  means 
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of  detecting  transition  location.  Cameras  mounted  in  a  wing  pod,  as  well  as  in  die 
cabin,  will  enable  DLR  and  Oxford  each  to  accomplish  transition  measurements  on 
both  sides  of  the  nacelle.  Drag  on  die  external  nacelle  will  be  assessed  by  DLR  using 
pressure  rakes  mounted  downstream  of  the  nacelle. 

The  instrumented  nacelles  are  to  be  manufactured  to  the  specifications  of 
Rolls-Royce,  Oxford  University,  and  DLR  by  Hurel-Dubois  (UK)  ltd.,  of  Burnley, 
England.  The  design  condition  for  the  nacelle  is  M  =  0.6  at  30,000  ft.,  althou^  drag 
and  transition  location  will  be  measured  diroughout  a  prescribed  fli^t  envelope. 

1.2  Transition 

The  transition  from  laminar  to  turbulent  flow  represents  one  of  the  most 
irr^rtant  and  perplexing  problems  in  fluid  mechanics.  REYNOLDS  (1883)  initiated 
the  study  of  transition  more  than  a  century  ago  with  his  classic  series  of  dye  filament 
experiments.  REYNOLDS  observed  that  a  diread  of  dye  flowing  down  a  pipe,  while 
initially  sharply  defined,  sutklenly  difiiised  into  the  stream  at  what  came  to  be  known 
as  the  location  of  ’transition’  from  lamiirar  to  turbulent  flow.  Furthermore,  he 
established  that  this  i^ienomenon  occuned  at  a  unique  numerical  value  of  die 
dimensionless  parameter,  U  dN  ~  the  transition  Reynolds  number. 

Since  that  time,  a  vast  amount  of  research  has  been  directed  at  better 
understanding,  predicting  and  controlling  die  phenomenon  of  transition  from  laminar 
to  turbulent  flow,  paiticulaiiy  within  boundary  layers.  Neverdieless,  a  complete 
understanding  of  boundary  layer  transition  remains  elusive. 
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1.2.1  A  Model  of  TroHsitioH 

There  currently  exists  no  goieralised  theory  which  would  find  its  foundation 
in  a  coiiq>rebensive  t^teciation  of  the  jAysical  mechanisms  responsible  for  transition 
and  extend  to  a  prediction  of  transition  onset  and  lengdi.  Nevertheless,  many  links 
in  the  chain  of  events  involved  wifo  transition  have  been  exposed.  In  particular,  the 
initial  stage  leading  to  transition,  the  ’linear’  portion  of  transition,  is  now  fairiy  well 
understood.  REYNOLDS  (1895)  postulated  that  transition  is  a  consequence  of  the 
instability  of  the  laminar  flow.  However,  it  was  the  theoretical  studies  of  PRANDTL 
(1922)  which  eventually  led  to  a  ’linear  stability  theory’  cq>able  of  predicting  a 
’critical’  Reynolds  number  at  whidi  smaU  disturbances  in  the  flow  (ToUmien- 
Schlichting  waves)  become  unstaUie.  The  results  of  linear  stability  theory  were  later 
confirmed  experimentally  by  SCHUBAUER  and  SKRAMSTAD  (1948).  SCHLICHTI- 
NG  (1979)  provides  a  detailed  account  of  die  history  and  development  of  linear 
stability  theory.  Here  linear  stability  theory  (or  ’classical  small  perturbation  analysis’) 
will  only  be  mentioned  briefly  as  a  conveniem  means  of  qipioaching  the  processes 
of  transition  and  the  effects  of  the  disturbance  environment  on  transition. 

Stability  analysis  is  based  on  the  assunqition  diat  laminar  flows  are  affected 
by  certain  small  disturbances  ^lich  may  originate  from  surfoce  inqierfections  or 
irregularities  in  the  external  flow.  The  theory  investigates  the  bdiaviour  of  these 
disturbances  as  they  are  superimposed  on  the  Navier-Stdres  equations  of  motion.  If 
die  disturbances  decay  with  time,  the  flow  is  considered  stable.  Disturbances  vdiidi 
amplify  with  time  are  unstaUe  and  may  cause  transition.  The  objea  is  to  predict  the 
value  of  the  critical  Reynolds  number  for  a  prescribed  flow. 
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The  fundamental  differential  equation  for  a  distuibanoe  is  the  well  known  Orr- 
Sununeifcld  (O-S)  equation  (die  derivation  of  which  may  be  found  in  texts  sudi  as 
SCHUamNG  (1979)): 

((/-c)(<|>"-a*<»)-l/"<t»  =  — +  .  (1.1) 

Here,  R(  represents  the  Reynolds  number  based  on  boundary  layer  thickness,  S,  a  is 
the  wave  number,  is  the  anqilitude  function  of  the  fluctuation,  and  c  is  a  complex 
function  of  which  the  real  portion  denotes  the  velocity  of  propagation  of  the  wave 
while  the  imaginaty  portion  determines  the  degree  of  danqiing  (or  amplification)  of 
the  disturbance.  Equation  1.1  is  obtained  by  writing  each  dependent  variable  in  the 
Navier-Stokes  equations  as  a  dependent  variable  plus  a  time  dependent  term.  The 
resulting  equations  are  linearised  and  the  mean  flow  equation  is  subtracted  out. 
Parallel  flow  is  assumed  and  the  equations  are  nondimensionalised  with  hee-stmam 
quantities.  Suitable  boundary  conditions  can  dien  be  inqiosed. 

By  omitting  die  viscous  terms  (RHS)  of  the  O-S  equation,  RAYLEIGH  (1913) 
was  able  to  derive  the  important  ’point  of  inflection  criterion’  which  states  that  the 
existence  of  a  point  of  inflection  in  the  velocity  profile,  U(y)  constitutes  a  necessary 
condition  for  the  occurrence  of  instability.  TOLLMEEN  (1936)  later  showed  diat  such 
an  inflection  also  constitutes  a  student  condition  for  die  anqilificadon  of  disturb¬ 
ances. 


The  solution  to  the  O-S  equation  for  a  prescribed  laminar  flow  is  a  curve  of 
neutral  stability  as  shown  in  Rg.  1.5.  The  poim  on  the  curve  at  whidi  die  Reynolds 
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nnmber  has  its  smallest  value  (tangent  parallel  to  the  veitical  axis)  is  of  particular 
interest  since  it  indicates  the  limit  of  stability.  Below  diat  critical  Reynolds  nnmber, 
all  individual  oscillations  decay;  for  Reynolds  numbers  >  Re^  distnibances  at 
at  least  some  frequencies  are  anq>lified.  It  is  wotfo  noting,  however,  foat  Re^  is  not 
the  same  as  the  Reynolds  number  at  transitkm  start.  Re,  (the  experimental  criticd 
Reynolds  number).  Re^  lepiesents  the  Reynolds  number  at  which  some  individual 
distuibances  will  begin  to  anqtlify.  As  will  be  discussed  below,  the  culmination  of 
the  anqtlification  process  is  turbulence.  The  point  at  which  isolated  patches  of 
turbulence  first  begin  to  appeal  is  designated  Re^  The  anqtlification  process,  however, 
takes  stMiie  time.  Therefore,  Re,  must  be  downstream  of 


The  relationshq>  between  Re^  and  Re,  is  not  yet  well  understood.  More 
specifically,  no  theoretical  model  exists  which  is  capable  of  explaining  die  means  by 
which  two-dimensional  Tollmien-Schlichting  (T-S)  waves  eiqierience  non-linear 
anqilification,  eventually  leading  to  turfoulenoe.  The  temporal  and  spatial  scales 
required  for  this  anqdificadon  may  be  influenced  by  a  number  of  environmental 
factors  (discussed  in  Sectkai  1.2.2.).  The  means  by  udiidi  die  boundary  layer 
responds  to  these  factors  is  little  understood.  Consequendy,  a  reliable  method  for 
prediction  of  transhifm  location  is  still  a  utopia. 

Studies  involving  this  ’mm^inear’  region  of  instability  began  soon  afier  the 
aforementioned  woik  of  SCHUBAUER  and  SKRAMSTAD  (1948).  From  observa¬ 
tions  of  the  flow  on  a  water  taUe,  EMMONS  (1951)  postulated  dutt  transition  is  not 
initialed  along  a  conthnioos  front  but,  radier,  involves  a  transition  region  in  which 
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'spots’  of  tuibulence  form  and  expand  as  they  move  downstream.  Fully  tuibulott  flow 
occurs  where  these  turbuloit  spots  merged  to  cover  die  entire  qian  of  the  surface. 
Passage  of  the  spots  over  points  on  die  surface  results  in  alternations  of  laminar  and 
turtnilent  flow.  The  alternations  may  be  quandtadvely  described  an  intermittency 
factor,  y,  whidi  represents  the  fraction  of  time  any  point  spends  in  turbulent  flow. 
Hence,  Y=  0  at  the  onset  of  transition  and  Y=  1  when  the  transition  to  tutbulem  flow 
is  complete.  SCHUBAUER  and  KLEBANOFF  (19SS)  provided  details  of  the 
geometry  of  turbulent  spots  and  their  growth  as  diey  move  downstream  by  using  an 
electrical  spark  to  introduce  turbulent  spots  into  a  laminar  boundary  layer  (Fig.  1.6). 
NARASIMHA  (19S7)  postulated  diat  the  spots  formed  at  a  preferred  stream-wise 
position  but  randomly  in  time  and  span-wise  location.  The  resulting  intermittency 
distribution  dirougfl  the  transition  region, 

^  (1*2) 

where 

r  =  ,  (IJ) 

A. 

fits  the  experimental  data  well.  Here,  A.  rqnesents  the  distance  between  the  location 
at  which  y  =  0.2S  and  the  location  at  which  y  =  0.7S,  and  x,  marks  the  beginning  of 
transition.  DHAWAN  and  NARASIMHA  (1958)  showed  diat  the  velocity  profile  and 
skin  frictitm  variations  widiin  the  transition  region  at  subsonic  speeds  could  be 
predicted  by  weig}iting  die  average  of  die  laminar  and  turbulent  values  by  die 
intermittency  factor.  The  validity  of  this  proposition  has  reoendy  been  supported  by 
die  results  of  DEY  and  NARASIMHA  (1990). 


SCHUBAUER  (1957)  and  KLEBANOFF  and  TEDSTROM  (1959)  began 
studies  into  the  manner  in  which  these  tnibolent  q>ots  arise  from  die  two-dimensioiial 
distuifoances  in  the  boundary  h^er.  These  studies  indicated  diat  diere  is  a  definite  and 
reproducible  progresskm  of  events  by  T-S  waves  evolve  into  turbulence.  Using 
a  vibrating  ribbon  technique  for  die  introduction  of  forced  oscillations  in  the  boundary 
layer,  KLEBANOFF  et  al.  (1962)  were  able  to  show  diat  even  in  nominally  two- 
dimensional  flow,  the  T-S  waves  that  develop  in  regions  of  instability  increasingly 
warp  or  bend  in  the  span-wise  ditectitm  as  they  move  downstream.  This  warping  is 
responsible  for  secondary  flows  to  and  from  the  surface  uhich  eventually  evolve  into 
turbulent  spots. 

Many  features  of  diis  evolution  from  T-S  waves  to  turbulent  spots  have 
become  widely  accepted  (see  Fig.  1.7).  STUART  (1965)  and  TANI  (1969)  provide 
general  reviews  of  the  transition  process.  It  nuiy  be  summarised  as  follows: 

(1)  Stable  laminar  flow  develops  a  region  of  instability  to  small  wavy  disturb¬ 
ances.  Provided  the  Reynolds  number  is  above  a  critical  level  (as  predicted  by  linear 
stability  theory),  T-S  waves  begin  to  anqdify. 

(2)  Three  dimensional  waves  and  vortices  form.  Minor  irregularities  in  die 
free  stream,  test  surface,  etc.  give  riae  to  a  rato  of  wave  gtowdi  which  varies  with 
span-wise  position.  Thus,  a  shown  in  Fig.  1.8,  an  initially  two-dimensional  wave 
takes  on  three-dimensional  fonn.  Vortex  filametts  widiin  the  boundary  layer  distort 
with  the  wave  front  The  vortex  filaments,  in  turn,  induce  secondary  flows  to  and 
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from  the  surface  which  help  to  intensify  die  waiping.  As  die  generated  diiee- 
dimensional  waves  with  voiticify  progress  downstream,  die  wave  front  distortitm 
increases.  The  disturbance  lines  dius  form  ’A*  or  ’hairpin’  vortices  (Fig.  1.9).  As 
these  grow,  the  vorticify  Unes  stretch  and  strengthen.  ’Peaks’  and  ’hollows’  of  high 
and  low  disturbance  intensify  form  (see  Fig.  1.10). 

(3)  Bursts  of  turbulence  appear  in  regions  of  very  high  vorticity.  The 
instantaneous  velocity  profile  develops  a  region  of  hi^  shear  (with  an  associated 
inflectitm)  at  span-wise  stations  cotiespcmding  to  a  peak  (Fig.  1.11).  As  the  point  of 
inflection  criterion  suggests,  these  inflectional  velocity  profiles  ate  associated  with 
higher  anqilification  rates.  The  vortices  act  as  a  mixing  mechanism,  introducing  low 
speed  flow  from  die  lower  layers  of  the  boundary  layer  to  the  upper  layers  and  vice 
versa.  Thus,  additional  tegicHis  of  inflectional  velocity  profiles  are  created,  dtereby 
adding  additional  sources  of  instability  and  disturbance  growth.  At  some  stage  high 
frequency  bursts  of  mrbulence  (embryonic  ’turbulent  spots’)  form  in  the  regions  of 
high  vorticity. 

(4)  Turbulent  spot  growth.  The  embryonic  turbulent  qiots  travd  downstream 
and  spread  span-wise  at  rates  as  observed  by  SCHDB AUER  and  KLEBANOFF  (1955) 
(refer  to  Fig.  1.6).  The  initiation  of  turbulent  ^x>ts  is  generally  considered  to  be  die 
’start’  of  transition.  At  some  distance  downstream,  these  qiots  coalesce  to  fr»m  a 
fully  turbulent  boundary  layer. 


The  details  of  transition  presented  above  are  primarily  derived  from 


’controlled’  experiments  (e.g.  diose  in  which  the  boundary  layer  is  disturbed  in  a 
controlled  manner)  such  as  those  performed  by  KLBBANOFF  et  al.  (1962)  employing 
a  vibrating  ribbon.  Nevertheless,  the  features  discussed  are  thought  to  be  common  to 
most  forms  of  ’natural’  transition  as  well,  althou^  the  events  pethiq>s  occm  in  a  less 
tegular  manner.  In  some  cases,  however,  it  seems  that  the  disturbance  enviromnent 
may  be  so  intrusive  upon  the  boundary  layer  as  to  ’bypass’  die  initial  T-S  wave  stage 
of  the  transition  process  and  introduce  directly  disturbances  of  large  enough  amplitude 
to  readily  generate  turbulent  spots.  For  exanq>le,  SUDER  and  O’BRIEN  (1988)  and 
SOHN  et  al.  (1989)  have  repotted  bypass  transition  due  to  hi^  levels  of  free  stream 
turbulence.  MORKOVIN  (1984)  has  also  addressed  the  issue  in  depth. 

1.2,2  Factors  Affecting  Transition 

Most  of  the  difficulty  in  understanding  and  predicting  transition  is  due  to  the 
large  number  of  independent  factcurs  whidi  affect  transition.  It  is  known,  for  exanqtle, 
that  transition  may  be  influenced  1^  free  stream  turbulence,  pressure  gradient,  heat 
transfer  between  the  surface  and  the  fluid,  surBue  roughness,  surface  curvature,  Mach 
number,  Reynolds  number,  angle  of  attadt,  sweep  angle,  acoustic  radiatimi,  surface 
vibration  and  injection  or  suction  of  fluid  at  the  wall.  STUART  (1965)  and  TANI 
(1969)  and,  more  recently,  RESHOTKO  (1976)  and  HERBERT  (1988)  have  produced 
excellent  general  reviews  of  transitkm  and  die  various  factcns  affecting  it*.  The  more 

*STUART  and  TANI  focus  on  low  speed  flows.  RESHOTKO  is  primarily 
concerned  with  heating  and  die  receptivity  of  die  boundary  li^r,  i.e.  the  means 
which  a  particular  forced  disturbance  enters  die  boundary  Iqw  and  die  nature  of  its 
signature  in  die  disturbance  flow.  HERBERT  is  inteiested  in  the  means  by  whidi 
small  two-dimensioital  distnrbaiices  in  flat  {date  boundary  layers  grow  into  the  violent 
diree-dimensional  disturbances  and  hi^-firequency  motions  that  diaracterise 
turbulence. 
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general  worics  of  SCHUCHTING  (1979),  YOUNG  (1989)  and  WHITE  (1991)  also 
discuss  the  various  factors  affecting  transititHi. 


The  effects  of  heating  and  pressure  gradient  on  transition  are  of  pritnary 
interest  in  the  present  study.  Curvature,  surface  rou^uiess  and  turbuleitce  are  of 
subsidiary  interest  and  will,  consequently,  be  discussed  briefly.  The  remaining  factors 
will  not  be  discussed  further. 


1.2,2.1  Curvature 

While  of  little  bearing  for  the  flat  plate  experiments  presented  in  diis  study, 
curvature  may  be  of  interest  for  experiments  performed  on  the  NLF  nacelle. 
LIEPMAN  (1943)  has  shown  that  Re,  is  virtually  unaffected  by  convex  curvature  when 
the  radius  of  curvature  is  large  in  relation  to  the  boundary  layer  thickness,  5  (Fig. 
1.12). 


1.2JI.2  Free-stream  Turbulence 

Qualitatively,  trartsition  is  hastened  by  an  increase  in  fiee-stream  turbulence. 
Most  studies  involved  with  turbulence  to  date  have  only  measured  fiee-stream 
turbulence  intensity. 


Tit  = 


+  v'*  +  w'*  , 


U 


1 

2 


(1.4) 


where  primes  denote  perturbed  quantities.  In  recent  years  there  has  been  a  growing 
awareness  of  the  importance  of  other  measures  of  turbulence  sudi  as  die  velocity  and 
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diiection  of  disturbance  propagation  as  well  as  tbe  frequency  and  wavelength  of  the 
disturbance.  MEIER  et  al.  (1987),  for  exaiiq>le,  have  shown  diat  the  <niset  of 
transition  correlates  better  with  the  spectral  density  of  the  free-stream  turbulence  in 
the  frequency  range  of  the  T-S  waves  than  the  total  intensity. 

The  effect  of  Tu  on  transition  locatirm  as  observed  by  DRYDEN  (1959)  is 
shown  in  Fig.  1.13.  There  appears  to  be  no  significant  effect  of  external  turbulence 
on  Re,  for  turbulence  levels  of  less  than  -  0.1  %.  With  further  increase  in  Tu,  diere 
is  a  continuous  decrease  in  Re,.  The  experimental  results  of  BLAIR  (1982)  and 
GOSTELOW  and  BLUNDEN  (1988)  show  a  similar  decrease  in  Re,  with  increasing 
Tu  above  a  turbulence  intensity  of  0.1%.  VAN  DRIEST  and  BLUMER  (1963)  and 
McDonald  and  FISH  (1973)  have  predicted  theoretically  the  same  trends.  HALL 
and  OroBINGS  (1972)  and  ABU-GHANNAM  and  SHAW  (1980)  have  provided 
further  experimental  data  atxl  empirical  predictions  accounting  for  fiee-stream 
turbulence  intensity  and  pressure  gradient.  EDWARDS  and  FURBER  (1956)  and 
JUNKAN  and  SEROVY  (1967)  have  found  no  significant  influence  of  7V<  on 
convective  heat  transfer  from  a  flat  plate,  apart  from  that  due  to  an  advancement  of 
transition  location. 

Ettqrirical  correlations  currently  provide  the  best  prediction  of  transitirm 
location  widi  variations  in  turbulence  intensity  and  pressure  gradient.  Figure  1.13,  for 
example,  shows  diat  the  semi-ertqmcal  prediction  formulated  by  MACK  (1977)  fits 
DRYDEN’s  data  quite  well  at  rtKxJeiate  tutbulence  levels.  At  very  low  totbulaice 
levels  (i.e.  Tu  <  0.1%),  however,  contenqrorary  onpirical  methods  frul  to  fit  die 


experimental  data  well.  Recently.  GOVINDARAJAN  and  NARASIMHA  (1991)  have 
postulated  that,  at  these  tuibulence  levels,  the  dominant  factor  influencing  the  onset 
of  transition  is  no  longer  fiee-stieam  turbulence,  but  residual  disturbances  (e.g.  noise, 
and  vibration)  specific  to  the  each  facility  and,  peifa^,  each  operating  condition. 
This  is  consistent  with  the  experimental  observations  of  DRYDEN  mentioned 
previously.  The  in^lication  is  that  die  mysterious  ’unit  Reynolds  number  effect’  may 
disqipear  provided  adequate  allowances  are  made  for  variations  in  these  sort  of 
parameten.  GOVINDARAJAN  and  NARASIMHA  propose  an  ’equivalent  turbulence 
intensity’,  q^,  to  represent  these  residual  facility-dependent  disturbances.  Their 
correlation  employing  presented  in  Fig.  1.14  against  a  range  of  experimental  data. 

1.2J1.3  Sutface  Roughness 

Surface  rou^iness  above  a  certain  scale  is  know  to  move  transition  forward. 
Two  critical  excrescence  heights,  denoted  k„i,j  and  k^,^,  are  of  interest.  Transition 
is  unaffected  by  excrescences  of  heights,  k,  where  k  <  k^,j.  If  k  >  k^,j,  increasing 
k  moves  transition  forward  until,  at  k  =  k^,j,  transition  occurs  just  af)  of  the  excres¬ 
cence.  The  values  of  k^,j  and  k^j  for  a  given  flow  depend  on  the  fluid  viscosity  and 
external  velocity  distribution  as  well  as  the  shiqie  and  location  of  individual 
excrescences. 

YOUNG  and  PATERSON  (1981)  provide  a  detailed  account  of  die  effects  of 
excrescences  on  traiuition  as  well  as  calculation  procedures.  The  only  theoretical 
study  known  to  dte  author  which  attempts  to  iiKlude  the  influence  of  surface 
roughness  on  transition  location  is  that  of  McDONALD  AND  FISH  (1973). 
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1.2^.4  Heating  and  Pressure  Gradient 

It  is  well  established  that  heat  transfer  to  or  from  a  surface  influences  boundary 
layer  transition.  LINKE  (1942)  observed  diat  heating  a  body  in  air  increased  its  drag, 
and  he  correctly  surmised  that  the  transition  location  had  been  moved  upstream.  The 
destabilising  influence  of  heating  in  air  is  chiefly  due  to  the  variation  of  viscosity  with 
temperature.  As  discussed  below,  hitter  viscosity  near  the  hot  surface  leads  to  an 
inflection  in  the  velocity  profiie  (Fig.  1.15a),  as  confirmed  experimentally  by 
LIEPMAN  and  FILA  (1947).  Cooling  leads  to  ’fuller’  velocity  profiles  and,  therefore, 
greater  stability  and  a  delay  of  transition  location  in  an  air  boundary  layer  (see 
MACK,  1975).  As  STRAZISAR  et  al.  (1977)  observed,  the  trends  ate  reversed  in 
water.  A  number  of  theoretical  investigations,  such  as  those  performed  by  WAZZAN 
et  al.  (1968),  NAYFFH  and  EL-HADY  (1980)  and  LEKOUDIS  (1980),  have  predicted 
the  same  trends.  WAZZAN  et  al.  have  noted  that  there  exists  a  critical  heating  rate 
for  water  flows.  Above  this  critical  heating  rate,  the  stabilising  influence  of  ’fuller’ 
velocity  profiles  is  overcome  by  the  destabilising  influence  of  reduced  viscous 
damping  and  a  further  increase  in  heating  reduces  the  Re^,. 

The  influence  of  pressure  gradient  on  transition  is  also  well  documented.  For 
example,  VAN  DRIEST  and  BLUMER  (1963),  HALL  and  GIBBINGS  (1972),  ABU- 
GHANNAM  and  SHAW  (1980)  and  BLAIR  (1982)  dealt  with  pressure  gradient  and 
free-stieam  turbulence  effects  on  the  transition  process.  A  favourable  pressure 
gradient  leads  to  delayed  transition  and  an  adverse  gradient  to  early  transition.  These 
effects,  like  those  due  to  heating,  may  be  explained  by  reference  to  velocity  profiles 
and  the  point  of  inflection  criterion  (see  Fig.  1.15b). 
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Neglecting  buoyancy  effects,  but  accounting  for  the  temperature  dependency 
of  the  viscosity  of  the  fluid,  n,  the  two-dimensional  boundary  layer  momentum 


equation  (Eqn.  13.6  in  SCHUCHTING  (1979))  may  be  written  as 

= .  (1.5) 

dx  dy)  dx  dy[  dy) 

where  x  is  distance  along  the  surface,  y  is  distance  normal  to  the  surface,  p  is  the 
mean  static  pressure  and  u  and  v  are  components  of  the  velocity  in  the  jc  and  y 
directions  respectively. 


Evaluated  at  the  wall,  u  =  v  =  0.  Hence, 


dx  9y(*^  0yj 


Expanding, 


where  w  denotes  values  evaluated  at  the  wall.  Rearranging  Eqn  1.7, 
{  d^u  _ _ ^  ^  ^ 


(1.6) 


(1.7) 


(1.8) 


where  the  LHS  of  Eqn  1.8  represents  the  curvature  of  the  velocity  profile,  «(y),  at  the 
wall.  If  the  wall  is  heated,  the  tenqrerature  of  the  fluid  at  the  wall  will  be  greater  than 
the  temperature  of  the  fluid  in  the  flee  stream.  The  viscosity  of  air  increases  with 
temperature  and,  therefore,  idft/dyX,  <  0  for  a  heated  surface  in  air.  SirKe  (du/dy).  > 
0  in  a  boundary  layer,  the  curvature  of  the  velocity  proflle  at  the  wall  is  positive  when 
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dp/dx  =  0.  The  cuivature  at  the  edge  of  the  boundary  layer,  however,  is  itegative.  It 
follows  that  a  point  of  inflection  exists  widiin  the  boundary  layer  velocity  profile. 


When  a  favourable  pressure  gradient  {-dptdx)  is  combined  with  surface  heating 
the  sign  of  the  velocity  profile’s  curvature  at  the  waU  depends  on  the  relative 
magnitudes  of  the  teims  on  the  RHS  of  Equation  1 .8.  The  boundary  layer  is  expected 
to  be  more  stabile  in  relation  to  the  boundary  layer  on  an  unheated  flat  plate  with  zero 


pressure  gradient  when 


Iflii)  (in)  <  ^1 


Conversely,  the  boundary  layer  is  relatively  less  stable  when 


To  the  author’s  knowledge,  there  ate  no  studies  in  the  literature  which  have 
quantitatively  examined  the  combined  effects  of  heat  transfer  and  fevourable  pressure 
gradient  on  transition  locatitm  in  an  air  boundary  layer.  CEBECI  and  SMITH  (1968) 
unsuccessfully  attempted  to  do  so  while  endeavouring  to  use  heating  as  a  boundary 
layer  tr^.  BLAIR  (1982)  examined  the  influence  of  pressure  gradient  on  transititm 
on  a  heated  surface  but  did  not  examine  die  influence  of  varying  heating  rates. 
LAUCHLE  and  QURNEY  (1984)  examined  the  effects  of  heating  on  an  underwater 
axisymmetric  body  with  a  neady-unifbrm  favourable  pressure  gradient  The 
temperature  distributirm  was  neatly  constant,  and  transition  location  was  plotted  versus 
heating  as  shown  in  Fig.  1.16.  However,  due  to  the  aforementioned  conqieting 


influences  of  viscosity  on  the  stability  of  a  water  boundary  layer,  quantitative 
agreement  on  the  movement  of  transition  with  heating  for  air  and  water  boundary 
layers  seems  unlikely.  The  results  presented  herein  represent  an  investigation  to 
examine  quantitatively  the  effects  of  heating  on  transition  location  with  both  zero  and 
favourable  pressure  gradients  in  die  incompressible  air  boundaiy  li^er. 

1.23  Techniques  for  Transition  Detection 

The  wholesale  change  in  the  characteristics  of  the  boundary  layer  flow  during 
transition  lends  itself  to  a  number  of  detection  techniques.  Some  of  the  most  common 
will  be  discussed  here.  A  study  of  the  velocity  distribution  in  the  boundary  layer,  for 
example,  offers  a  conveniem  means  of  detecting  transition.  Fig.  1.17  shows  that 
transition  involves  a  significant  change  in  the  shape  of  die  mean  velocity  distribution 
curve  through  the  boundary  layer.  It  is  also  appazetit  from  this  figure  that  the 
boundary  layer  suddenly  thickens  in  the  region  of  transition.  The  flat  plate  boundary 
layer,  for  exaitq)le,  changes  from  a  laminar  growth  rate  pn^rtional  to  to  a 
turbulent  growth  proportional  to  x^.  The  shape  factor,  H  =  5/5}  also  decreases  from 
>>2.6  in  die  laminar  region  to  »  1.4  in  die  turbulent  region  (Fig.  1.18).  Here,  Sj  is 
the  displacement  thickness  and  8^  the  momentum  thickness.  Hot  wires,  hot  films  and 
Pitot  tubes  offer  convenient  means  of  measuring  the  boundary  layer  velocity 
distribution. 

The  higher  velocity  gradient  at  the  wall  associated  with  turbulent  boundary 
layers  is  responsible  for  increased  shearing  stress,  t  according  to  the  foimula. 
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The  characteristic  rise  in  the  coefficient  of  skin  friction  through  the  transitional  region 
is  amenable  to  detection  by  a  number  of  techniques  including  surface  shear  stress 
gauges  and  shear  sensitive  liquid  crystals  (on  which  more  will  follow). 

In  the  transition  region,  turbulent  mixing  leads  to  a  rise  in  heat  transfer  and 
mass  transport  which  is  analogous  to  the  rise  in  skin  friction.  This  increased  mass 
transport  is  utilised  by  sublimttting  chemicals  to  indicate  transition  location.  The 
chemical  sublimes  first  at  the  point  of  highest  mass  transport  (transition  end)  and  last 
at  the  point  of  lowest  mass  transport  (transiticMi  start).  Thermocouples  or  terrqjerature 
sensitive  liquid  crystals  (on  whidi  mote  will  be  said  later)  may  be  enqrloyed  to  detect 
temperature  variations  on  a  uniformly  heated  surface  in  unheated  flow,  or  an  unheated 
surface  in  heated  flow.  DEOM  (1989),  CROWDER  (1990)  and  HORSTMAN  et  al. 
(1990)  have  recently  reported  on  thermo-imaging  as  a  non-intrusive  means  of  mapping 
surface  tettqwratute.  Given  a  knowledge  of  die  flow  characteristics  and  die 
tenqierature  distribution,  transition  start  and  end  may  be  determined. 

As  noted  in  Section  1.2.1,  intermittency  also  provides  a  means  of  measuring 
the  transition  location.  Fig.  1.19  shows  hot  film  traces  at  various  locations  dirough 
an  intermittent  (transitional)  boundary  layer.  Intermittency  increases  from  a  value  of 
Y  ">  0.0  in  the  laminar  region  to  y  »  1.0  at  the  end  of  transition.  An  ’intermittency 
meter’  such  as  that  cwistructed  by  KIM  et  al.  (1989),  may  be  used  to  translate  hot 
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film  intennittency  traces  into  quantitative  values  of  y. 

Intermittency  is  generally  recognised  as  providing  die  best  indication  of  die 
location  of  the  transitional  region  since,  jdiysically,  the  location  at  which  tuibulent 
spots  first  appear  is  represoited  by  y  =  0.0,  and  1-0  coitespcmds  to  the  location  at 
which  the  sfmts  coalesce  to  form  a  fully  turbulent  boundary  layer.  OWENS  (1970) 
and  BLAIR  (1982)  have  noted,  however,  that  other  methods  of  detecting  transition 
may  not  provide  the  same  transition  start  and  md  locatioos  as  those  indicated  by 
intermittency.  This  is  to  be  expected  if,  as  NARASIMHA  (1958)  suggested,  mean 
flow  characteristics  in  the  transitional  region  may  be  considered  to  be  a  combination 
of  fiilly  laminar  and  fully  turbulent  flows,  each  pn^rtitmally  weighted  by  the 
intermittency  faaor.  At  low  intennittency  factors,  the  turbulent  contribution  to  the 
mean  flow  characteristics  will  be  small.  Skin  biction,  for  example,  will  cmitinue  to 
fall  for  a  short  period  after  Y  >  0.  Hence,  die  transitimi  location  as  indicated  by 
minimum  t  will  be  downstream  of  the  location  at  which  the  value  of  y  has  begun 
increasing. 

Nevertheless,  die  difference  between  transition  locations  as  derived  from  heat 
transfer  distributions,  for  example,  and  those  determined  from  intermittency  is  small. 
GOVINDARAJAN  and  NARASIMHA  (1991)  tentatively  suggest 

Re',  -  0.9  Re,  ,  (1-12) 

where  Re’,  is  the  Reynolds  number  at  which  turbtdent  qiots  first  appear  and  Re,  is  that 
indicated  from  the  heat  transfer  or  skin  friction  distribution. 
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13  Liquid  Crystals 

Liquid  aystals  (LC's)  are  a  higjbly  anisotnq>ic  state  of  matter  wfaidi  exists 
between  the  solid  and  liquid  {biases  of  some  conqwunds.  Aldunigh  diey  qipear  as 
liquids,  they  exhibit  optical  prqreities  characteristic  of  a  crystalline  (solid)  state.  In 
particular,  LC’s  scatter  li|^  very  selectively. 

This  behaviour  is  trtu^ed  to  their  molecular  arrangement.  Each  LC  has  a  helical 
stmcture  whose  characteristic  pitch  length  fortuitously  coincides  with  the  wavelength 
of  the  visible  spectrum.  When  the  heUx  is  subjected  to  certain  physical  influences  (the 
two  primaiy  influences  being  temperature  and  siuface  shear  stress),  the  helix  pitch 
changes  and  the  wavelength  of  reflected  light  is  altered  accordingly.  Since  the  LC’s 
fundamental  chemical  structure  is  unaffected  by  dtese  changes,  its  responses  are 
theoretically  rapid,  continuous  and  reversible  indefinitely. 

Measurements  of  a  particular  physical  mechanism  are  made  by  ensuring  foat 
die  LC’s  sensitivity  to  odier  stimuli  is  either  reduced  to  a  ne^igible  level  or  accounted 
for.  For  diese  tests,  encipsulated  chiral  nematic  LC’s  were  used  for  temperature 
indication.  Hie  gel-like  ctpsule  enveloping  eadi  LC  consists  of  polyvinyl  alcohol  and 
is  apfuoximately  10  pm  in  diameter.  In  this  form,  the  LC’s  sensitivity  to  shear  is 
ne^igiMe  and  surface  temperatures  can  be  measured  to  within  0.2  °C.  IRELAND  and 
JONES  (1987)  rqxHt  re^ioiise  times  of  -  3  ms  for  these  eireqisulated  LC’s.  A  neat 
(unencqisulated)  cholesteric  liquid  crystal  whidi  has  a  temperature  insensitive 
selective  colour  reflection  was  used  for  shear  indication. 
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Over  the  past  20  years  LC  coatings  have  found  various  uses  in  aerodynamic 
and  aerothermal  studies.  JONES  (1991)  discusses  the  various  properties  of  LC’s 
which  may  be  useful  to  research  engineers.  Only  those  liquid  crystal  tediniques 
which  are  relevant  to  the  present  experiments  are  addressed  here.  OKER  and  MERTE 
(1976),  HIPPENSTEELE  et  al.  (1983),  BAUGHN  et  al.  (1985,  1986),  and  FISHER 
and  EIBECK  (1990)  have  used  umperatwre  sensitive  LC  coatings  to  measure  die 
surface  temperature  in  steady-state  convective  heat  transfer  experiments.  IRELAND 
and  JONES  (1985,1986),  JONES  and  HIPPENSTEELE  (1988),  BYERLEY  et  al. 
(1988)  and  WANG  et  al.  (1990)  have  utilised  tenqieratute  sensitive  LC’s  to  measure 
convective  heat  transfer  using  a  transient  tedinique.  In  the  steady-state  tests,  a 
knowledge  of  the  power  being  supplied  to  a  heated  strip  and  the  surface  and  free- 
stream  tenqieratures  enables  a  heat  transfer  coefficient  to  be  determined.  In  transient 
tests,  the  heat  transfer  coefficiem  is  found  from  the  fiee-stream  temperature  and  the 
surface  tenqierature  history  as  indicated  by  the  LC.  For  tenqteiatute  indication, 
cholesteric  liquid  crystals  which  have  a  narrow  colour  play  tenqieratuie  band  (typically 
1  ”0  are  chosen.  The  tenqieiature  level  at  whkh  this  colour  play  occurs  can  be 
adjusted  over  a  wide  range  of  temperatures  by  modifying  the  composition.  A  direct 
conqMuison  of  the  steady-state  and  transient  techniqiies  by  BAUGHN  et  al.  (1989) 
concluded  that  the  fonner  was  sinqder  and  provided  lower  uncertainties,  while  die 
latter  was  better  suited  to  conqdex  geometries.  The  present  woik  is  intended  to  assess 
the  suitability  of  the  steady-state  teciiniqoe  qiecificaily  for  transitkm  detection. 

The  use  of  shear  sensitive  liquid  ciystals  in  aerodynamic  measuremaits  has 
been  a  tnoae  recent  knovatkin.  Two  difierent  processes  can  be  onployed  to  indicate 


shear  level  with  this  form  of  LC.  In  the  first  process,  the  LC  layer  starts  widi  a 
Grandjean  texture,  displaying  colour,  and  shearing  produces  a  colour  change.  This 
process  is  reversible  as  the  material  will  assume  its  original  colour  when  shearing 
ceases.  GALL  and  HOLMES  (1986),  HOLMES  and  OBARA  (1987),  and  JOHNSON 
(1988)  have  used  the  neat  cholesteric  LC  coatings  to  indicate  shear  stress  by  tiiis 
reversible  colour  play  technique  in  flight  tests.  REDA  (1991)  utilised  shear  sensitive 
LC’s  to  study  unsteady  fluid  physics  on  an  oscillating  airfoil,  reporting  a  response 
time  of  <  30  ms.  BONNET  et  al.  (1989)  note  that  the  reversible  technique  gives  a 
qualitative  measure  of  relative  shear  stress  but  is  not  useful  for  quantitative 
measurements  due  to  the  strcmg  dependence  of  illumination  and  viewing  angles  on  the 
colours  seen. 

In  the  second  process,  that  used  in  die  present  experiments,  the  material  is 
applied  to  the  surface  atxl  heated  in  order  to  drive  it  into  the  isotropic  phase.  The 
model  and  LC  film  are  then  allowed  to  cool,  the  LC  adopting  the  colourless  focal 
conic  texture.  If  this  is  then  sheared,  die  movement  of  the  film  aligns  the  cholesteric 
helices  and  the  layer  diqiiays  die  colourful  Grandjean  texture  (see  Bonnet  et  al. 
(1989)).  This  texture  diange  from  focal  conic  to  Grandjean  is  irreversible.  Such  a 
texture  change  mediod,  relying  on  the  time  to  change  from  a  colourless  to  a  colourful 
texture,  is  not  dependent  on  viewing  angle  or  illumination  and,  therefore,  may  be 
calibrated  to  give  a  quantitative  measure  of  shear  stress.  In  these  experiments,  the 
irreversible  texture  change  technique  has  been  used  for  the  first  time  to  obtain 
quamitative  measurements  of  surface  shear  stress  in  a  transitional  boundary  layer. 
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Liquid  crystals  o^r  a  number  of  potential  advantages  over  odier  transidtHi 
detection  techniques.  Liquid  crystals  are  simple  and  inexpensive  to  use.  Unlike  Pitot 
tubes,  hot  wires,  and  hot  films,  they  offer  fiiU  surface  m^ing  of  teiiq)erature  and 
shear  stress.  Reversible  LC  techniques  also  overcome  some  of  the  limitations  of 
sublimating  chemicals  by  providing  continuous  indications  of  transition  rather  than 
only  one  transition  pattern  per  flight.  Low  ambient  temperatures  confine  sublimating 
chemicals  to  in-flight  tests  below  20000  fi.,  while  HOLMES  et  al.  (1986)  have  used 
shear  LC’s  up  to  50000  ft.  Thermo-imaging  is  relatively  expensive  and  restricted  by 
surface  material,  surftu%  emissivity  and  conduction  effects  due  to  non-uniform 
substrates. 

Liquid  crystals  also  have  a  number  of  limitations.  The  heating  associated  with 
temperature  sensitive  LC  techniques  is  known  to  influence  the  transition  location.  The 
steady-state  heat  transfer  technique  is  unsuited  for  use  <hi  complex  geometries,  may 
be  time  consuming,  and  has  an  element  of  subjectivity,  although  AKINO  et  al.  (1989) 
have  developed  a  promising  means  of  alleviating  the  latter  two  difficulties  by 
employing  band-pass  colour  filters  and  a  blade  and  white  video  camera  in  conjunction 
with  the  LC’s.  Damage  to  the  LC  layer  due  to  ultraviolet  light  and  moisture  has 
previously  presented  an  obstacle  to  flight  tests  involving  temperature  sensitive  LC’s. 
WALTON  (1S)90),  however,  has  recently  ratted  progress  on  an  ultraviolet  blocking 
agent  as  well  as  a  waterproof  coating  for  temperature  sensitive  LC’s.  As  for  shear 
sensitive  LC’s,  die  reversible  colour  play  technique  offers  only  a  qualitative  indication 
of  surface  shear  stress.  While  the  irreversible  texture  change  technique  provides 
quantitative  measures  of  shear  stress,  the  crystal  layer  must  be  prqiared  anew  for  each 
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test  condition.  Except  for  the  irreversible  shear  sensitive  LC  technique,  LC’s  offer 
little  directional  flow  visualisation. 
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Chapter  2:  Experimental  Techniques 


2.1  Test  Facilities 

2.1.1  Tunnel 

All  CTpftrimfjits  were  carried  out  in  a  low  speed  wind  tunnel  (0  to  40  in/s) 
which  drew  in  air  at  ambient  pressure  and  tenqperature  (Fig.  2.1).  The  test  section  of 
the  tuntrel  measured  500  mm  wide  (nominally)  by  300  mm  in  height  For  each  metre 
advanced  downstream,  the  test  section  widened  by  qjproximately  20  mm  to  allow  for 
boundary  layer  grow*  along  the  tunnel  walls.  One  wall  of  the  test  seaion  consisted 
of  perspex  and  was,  therefore,  particularly  suited  for  viewiitg  liquid  crystals. 
Removable  secti<Mis  in  the  tunnel  roof  were  also  replaced  wi*  12  mm  perspex  in 
order  to  better  accommodate  viewing  and  i*otogrq.hing  die  test  surface.  H^’es  were 
drilled  in  diese  replacement  paireb  to  facilitate  various  probes  and  securing 
mechanisms  for  a  pressure  gradient  fairing. 

The  tunnel  was  not  paiticulariy  ’clean’.  Smoke  and  shear-sensitive  liquid 
crystal  flow  visualisation  teduiiques  suggested  diat  stwne  cross-flow  was  present  on 
a  flat  plfltf  insetted  in  flic  test  section.  Fiee-stream  turbulence  measurements  using 
a  constant  temperature  hot  wire  indicated  that  firee-stream  turbulence  intensity  was 
approximately  0.2  to  0.5  %.  Electrical  interference  generated  by  the  tunnel  motor 
precluded  more  accurate  measurements  of  flee-stieam  turbulence. 
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2.1  J  FlatPlaU 


A  symmetrical  flat  plate  served  as  the  measuremmt  surface.  It  was  machined 
firom  12  mm  perspex  and  measured  0.9  metres  in  length.  The  plate  was  secured  to 
two  aluminium  channeb  which  were  bolted  to  the  tunnel  floor.  The  q^aratus  was 
designed  to  produce  symmetrical  flow  on  either  side  of  the  plate.  T\innel  conditions 
could  be  adjusted  to  allow  natural  transition  to  occur  over  the  length  of  the  model. 
In  order  to  avoid  laminar  flow  separation,  a  leading  edge  involute  was  chosen  on  the 
basis  of  experimental  evidence  compiled  on  flat  plates  under  similar  flow  conditions 
by  MOSS  (1991)  and  DAVIS  (1980).  The  tq>er  ftran  the  leading  edge  expanded  to 
the  thickness  of  the  plate  at  SO  mm  downstream.  The  span  of  the  plate  increased 
slighdy  in  accordance  with  the  downstream  expansion  of  die  tunnel  width.  Machining 
of  the  plate’s  leading  edge  was  performed  on  an  NC  milling  machine.  A  seal  was 
formed  betwem  the  edges  of  the  plate  and  the  tunnel  wall  using  a  foam  sealing  sti^. 
Three  such  perspex  plates  were  made:  the  first  on  which  to  adhere  an  aluminium 
heating  strip,  a  second  with  pressure  tappings  of  0.5  mm  diameter  placed  at  ISO  mm 
intervals  along  its  centre-line  in  order  to  assess  pressure  gradients  along  the  surface 
of  the  perspex  plate,  and  a  third  in  which  to  mount  a  caitxm-fibre  hearar.  A  fourth 
plate  with  a  sligbtly  difierent  leading  edge  was  used  for  the  shear  sensitive  LC  tests. 

2.U  Pressure  Pairiug 

Favourable  pressure  gradkiits  tvere  qtfrikd  to  the  plate  using  an  adjustable  flat 
fairing  above  the  measurement  plate  as  shown  in  Figure  2.2.  Ibe  fairing  craisisted  of 
two  sections  (the  aft-  and  fore-fairings)  and  was  designed  to  produce  a  unifoim 
pressure  gradient  on  die  test  surface.  Bodi  sections  of  the  fairing  ass«nbly  were 
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constructed  of  12  mm  perspex  secured  to  the  tuimel  ceiling  and  hinged  at  a 
downstream  location  corresponding  to  the  leading  edge  of  the  measurement  surface. 
Thus,  the  aft-fairing  could  be  lowered  while  the  fore-fairing  remained  stationary.  The 
aft-fairing  extended  downstream  slightly  past  the  trailing  edge  of  die  measurement 
surface.  The  leading  edge  of  the  fore-fairing  tapered  into  the  tunnel  contraction.  A 
pressure  upping  in  the  fore-fairing  provided  a  means  of  measuring  oncoming  velocity. 
Curved  sections  attached  to  the  upper  surface  of  die  aft-fairing  protruded  through  the 
tunnel  roof  and  were  clamped  to  its  exterior,  locking  the  aft  fairing  in  position. 
Maridngs  on  the  curved  sections  facilitated  accurate  positioning  of  the  fairing.  The 
edges  of  the  fairing  were  sealed  with  foam  sealing  strip. 

2.2  Temperature  Sensitive  liquid  Crystal  Experimetds 

In  the  steady  state  convective  heat  transfer  experiments,  ohmic  heating  was 
used  to  supply  a  uniform  surftice  heat  flux  to  the  flat  plate.  Liquid  crystals  on  the 
surface  of  the  plate  were  then  used  to  measure  the  surface  tenperature  and,  given  a 
knowledge  of  flow  conditions  and  the  power  being  input  to  the  plate,  the  coefficient 
of  heat  transfer,  h,  was  determined. 

Transition  start  and  end  were  identified  by  the  characteristic  rise  in  mean  heat 
transfer  coefficient  in  the  transition  regitm  (Pig.  2.3a).  As  the  surftice  of  die  plate  was 
heated,  the  LC  would  display  contours  of  constant  temperature.  Qualitatively, 
transition  start  was  found  by  the  first  iqpearance  of  liquid  crystal  (LC)  colour  (lowest 
h  and  highest  surface  temperature)  as  diown  in  Fig.  2.3b.  As  die  current  flowing  to 
the  heater  was  increased,  the  isotherm  moved  upstream  in  the  laminar  region  and 


downstream  in  the  transitional  region.  An  isotherm  also  moved  upstream  through  the 
turbulent  region  until  it  eventually  merged  with  the  isotiterm  moving  downstream 
through  the  transitional  region  (Fig  2.3c).  This  merger  was  taken  to  indicate  the  end 
of  the  transitional  region  (highest  h  and  lowest  surface  temperature). 

As  noted  in  Section  1.2.3,  this  technique  does  not  necessarily  give  the  same 
transition  locations  as  those  determined  using  boundary  layer  shape  factors  or 
intermittency  techniques.  However,  this  method  of  identifying  transition  does  give  a 
consistent  set  of  data  on  transition  with  heat  transfer  and  pressure  gradient. 

For  these  experiments,  two  types  of  heating  strip  were  used.  An  aluminium 
heater  was  employed  to  validate  the  measurement  techiuque  and  establish  the  data 
reduction  procedures.  A  carbon  heater  was  then  used  to  simulate  that  which  is  to  be 
installed  on  the  NLF  nacelle. 

2.2.1  Aluminium  Heater 

2.2.1.1  Experimental  Apparatus 

A  127  mm  wide  heater  strip  was  adhered  along  the  centre-line  of  the  flat  plate 
as  shown  in  Fig.  2.4.  The  heater  str^  extended  from  the  downstream  end  of  the  plate 
around  the  leading  edge  and  back  to  the  rear  on  the  other  side  of  the  plate.  Electrical 
connections  were  made  at  each  end  of  the  strq}  at  the  rear  of  the  plate.  A  current  was 
passed  through  the  strip,  which  had  uniform  resistance,  to  supply  a  uniform  surface 
heat  flux  to  both  sides  of  the  plate. 
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The  heater,  of  50  |iin  overall  thickness,  consisted  of  a  polymer  film  metallised 
with  aluminium  and  was  su{^lied  by  Ooodfellow  ltd.  of  Cambridge,  England. 
Resistance  measurements  using  a  four-wiie  technique  showed  that  the  film  had  a 
resistance  of  0.46  fVsquare  at  18.S°C.  The  film’s  tenq)erature  coefficient  of 
resistance,  a,  was  measured  to  be  8.4  x  lO**  CiTC  by  placing  a  sample  of  the  film  in 
a  bath  of  de-ionised  water  and  measuring  its  resistance  over  a  40*^2  temperature  range. 
This  a  was  estimated  to  lead  to  maximum  variations  in  surface  heat  flux  of  <  0.26% 
for  the  transitional  boundaiy  layers  of  the  present  experiments.  OKER  and  MERTE 
(1981)  have  reported  that  the  resistance  of  thin  films  as  dqtosited  generally  changes 
with  age  as  the  strains  within  the  film  are  relieved.  This  annealing  process  is 
accelerated  by  heating  die  thin  film.  However,  no  measundile  change  in  resistance 
was  recorded  for  the  aluminium  film. 

A  qualitative  measure  of  the  unifoimity  of  the  heater’s  resistance  was  made  by 
passing  a  current  through  the  heater  under  zero  flow  and  diecking  for  unifoimity  in 
the  colour  play  of  the  LC’s  on  die  surface  of  the  heater.  Such  a  test  indicated  that  the 
resistance  was  quite  unifonn,  provided  the  section  of  aluminium  was  selected  with 
care.  Due  to  the  manufacturing  process,  resistance  was  found  to  decrease  near  one 
edge  of  the  roll  and  this  section  was,  therefore,  avoided.  No  quantitative  measurement 
was  made  of  the  unifoimity  of  the  heater’s  resistance.  BAUOHN  et  al.  (1985)  have 
reported,  however,  that  the  non-unifoimitMS  in  resistance  of  polyester  metallised  with 
vqiour-deposited  gold  were  found  to  be  less  than  5%  and  introduced  errors  in 
measurements  of  local  heat  transfer  coefficient  of  2-4%.  Using  similar  mateiiak, 
HIPPQ4STEELE  et  al.  (1983)  have  estimated  the  non-uniformities  in  resistance  to 
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cause  maximum  eirois  in  local  heat  flux  of  ±  8%,  although  this  estimate  could  be 
significantly  reduced  by  ’careful  selection’  of  the  heater  sheet. 

The  metallised  film  was  attached,  polymer  side  up,  to  the  plate  with  impact 
adhesive.  This  method  proved  simple,  produced  few  bubbles  between  the  substrate 
and  the  film,  and  made  for  an  adequate  bond  while  still  allowing  for  re-positioning 
when  necessary.  Care  was  taken  to  avoid  damaging  the  metal  side  of  the  film  as  any 
imperfections  quickly  caused  die  heater  to  fail  when  a  current  was  spplied. 

Good  electrical  contact  with  the  heater  was  found  to  be  inperative  as  points 
of  high  contact  resistance  quickly  burned  throu^  the  heater.  Adequate  contact  was 
achieved  by  using  a  metal  clanp  on  either  side  of  die  plate.  A  thin  strip  of  rubber 
was  sandwiched  between  the  polymer  side  of  the  heater  and  a  stiff  metal  strip. 
Similaily,  a  copper  strp  was  placed  between  the  metal  side  of  the  heater  and  the 
perspex  substrate.  Screws  then  ran  through  the  entire  assembly  and  were  threaded  into 
the  perspex.  By  tightening  the  screws,  the  aluminium  was  forced  securely  against  die 
copper  strips.  The  rubber  piece  betweoi  the  aluminium  strip  and  die  heater  made  for 
a  reduced  and  more  uniform  contact  resistance  between  die  tduminium  heater  and 
co(per  strip.  Wire  leads  soldered  onto  the  copper  strips  carried  electrical  currem  to 
the  assembly.  Figure  2.S  illustrates  dw  arrangement  of  the  electrical  contacts. 

A  very  thin  (12S  pm)  type  T  thermocouple,  attached  to  the  surface  of  the 
heater  film  near  die  rear  of  the  plate  (see  Fig.  2.4),  was  used  to  determine  the 
temperature  at  which  the  LC  changed  colour  (i.e.  to  calibrate  the  crystal).  The  LC 


displays  a  series  of  colours  corresponding  to  various  temperatures.  However,  only  the 
yellow  isotheim  was  calibrated  as  it  was  most  easily  identifiable  and  had  die 
narrowest  temperature  band.  The  yellow  isotherm  was  estimated  to  define  a 
temperature  to  within  0.2  °C.  Hie  upper  surface  of  the  heater  film  was  sprayed  with 
a  thin  coating  (approximately  S  pm)  of  a  matt  black  paint  to  provide  a  dark  back¬ 
ground  against  which  the  coloured  light  reflected  by  the  IX  layer  was  easily  viewed. 
The  painted  surface  was  smoothed  with  1200  grade  emery  pqier.  It  was  imperative 
that  the  painted  surface  be  allowed  to  dry  conqiletely  before  applying  the  LC  as  the 
trichloroethatie  solvent  contained  within  the  paint  is  damaging  to  encapsulated  LC. 

Thermochromic  enctqisalated  liquid  crystals  were  suf^lied  by  B.D.H.  ltd.  of 
Poole,  England.  The  LC’s  are  in  liquid  form  and  are  combined  with  Mowilidi  binder 
in  a  mixture  of  5  parts  LC  to  one  part  binder.  A  quantity  of  deminenilised  water 
equal  to  the  quantity  of  LC  and  binder  was  also  added  to  make  the  mixture  thinner 
and  easier  to  air  brush  onto  the  surface  in  thin,  uniform  coats.  Filtering  the  LC 
through  a  thin  mesh  screen  helped  prevent  die  air  gun  from  developing  blockages. 

Three  types  of  LC,  each  with  a  different  colour  play  temperature  range,  were 
pre-mixed  and  sprayed  onto  the  painted  surface  in  thin  layers.  This  resulted  in  a 
maximum  of  three  isotherms  being  produced  for  a  given  heat  flux  from  the  heater 
strqi.  Individual  layers  were  heated  with  a  hot  air  gun  to  evqiorate  moistuie  and 
access  the  quality  of  the  colour  play  (colour  intensity  increasing  widi  diickness  up  to 
a  limiting  value).  A  protective  layer  of  binder  was  spnycd  over  die  LC  layer.  This 
protective  layer  may  be  sanded  in  order  to  reduce  surface  roughness.  The  LC  and 
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binder  layers  were  q^roximately  20pm  thick.  Maximum  variations  in  the  surface 
finish  were  estimated  by  WALTON  (1990)  to  be  <  6  pm.  Calculations  based  on  a 
review  of  YOUNG  and  PATTERSON  (1981)  suggest  diat  diis  is  well  below  die  height 
at  which  a  two-dimensional  excrescence  would  disturb  the  boundary  layer.  The 
packaging  arrangement  of  the  heater  and  liquid  crystals  is  shown  in  Fig.  2.6. 

Power  was  supplied  to  the  heater  from  a  11  anq>,  240  volt  fiJC  variac 
transformer.  The  leads  to  the  heater  were  taken  through  the  roof  of  the  tunnel 
downstream  of  the  plate.  Current  was  measured  with  a  Salford  Selectest  model  ’Super 
SO’  aiudog  ammeter  which  was  connected  in  series  with  the  heater  and  power  supply. 
Current  was  used  to  calculate  die  power  being  dissqiated  by  the  heater.  These 
calculations  were  confirmed  using  voltage  measurements  frmn  a  Thuriby  model  1503- 
HA  digital  multimeter  which  was  connected  in  pataUel  with  the  heater.  The  electrical 
configuration  is  shown  schematically  in  Fig.  2.7.  Estimates  of  die  accuracy  of  all 
instruments  ate  listed  in  the  uncertainty  analysis,  ^rpendix  B. 

Mean  total  and  static  pressures  were  measiued  using  a  Pitot-static  tube  placed 
at  the  coitte-line  of  the  tunnel  at  mid  height  qipiaxiniaiely  10  mm  downstream  of  die 
plate.  The  Pitot-static  tube  was  connected  to  a  13  station  manometer  bank.  When  die 
plate  used  to  measure  surfme  pressure  gradioit  was  installed,  no  electrical  connections 
to  the  plate  were  necessary.  However,  five  additional  static  pressure  lines  ran  from 
the  static  ports  along  die  undoside  of  the  flat  plate.  The  lines  exited  throu^  the  roof 
at  the  rear  of  the  [date  and  were  connected  to  the  manometer  bank. 
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A  total  tempeiatoie  probe  insetted  through  die  tunnel  roof  alongside  the  ntot- 
stadc  tube  assessed  die  mean  total  tempeiatuie.  It  ccmsisted  of  a  type  K  diennocouide 
attached  to  a  ceramic  rod  and  housed  in  an  aspirated  aluminium  tube  of  l.S  mm 
diameter.  A  90°  elbow  held  die  thermocouple  normal  to  the  fiee  stream.  Torqierature 
measurements  were  displayed  digitaUy  <m  a  Comaik  mottel  6900  digital  anqdi- 
fierAhermometer. 

2Jt.l.2  Experimental  Program 

2.2.1.2.1  Flat  Plate,  Zero  Pressure  Gradient 

The  objectives  of  the  tests  using  ten^rature  sensitive  liquid  crystals  on  a  flat 
plate  at  aero  incidence  were  duee-fold:  (1)  to  ensure  that  natural  transition  would 
occur  on  the  measurement  surface,  (2)  to  verify  that  transition  could  be  detected  and 
analyzed,  and  develop  procedures  to  do  so,  and  (3)  examine  the  effects  of  heating  on 
transition  location  under  a  zero  pressure  gradient. 

Pressure  tappings  placed  at  ISO  nun  intervals  indicated  dut,  due  to  leakage  in 
the  tunnel  doors,  an  adverse  pressure  gradient  existed  along  die  plate  in  die  stream- 
wise  direction.  This  was  corrected  by  sealing  the  joints  of  the  tunnel  widi  foam  tiqie 
and  clanging  the  tunnel  at  qiproximately  die  centre  of  the  test  section  with  a  C- 
clanqi.  At  zero  incidence,  the  test  surface  then  had  a  zero  pressure  gradiem. 

Static  pressure  measuiements  on  die  test  surface,  when  conqiared  with  those 
of  the  Pitot-static  tube,  indicated  that  the  fiee-stream  flow  qieed  increased  slighdy  in 
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the  region  of  the  flat  plate  (see  Fig.  2.8).  Therefore,  measurements  of  actual  fiee- 
stream  velocities  in  the  region  of  the  plate  were  obtained  by  conecting  the  velocity 
measurements  provided  by  tite  Pitot-static  tube  accoiding  to  the  foimula, 

«  =  1.017  +0.39  (2.1) 

where  u  is  the  velocity  over  the  test  surface  and  u,  is  that  indicated  by  the  Pitot-static 
tube. 


An  AC  current  was  passed  through  the  heater  strip.  The  width  and  resistance 
of  the  strip  being  uniform,  die  current  produced  a  uniform  ohmic  heat  flux  along  the 
length  of  the  strq).  The  system  was  left  to  reach  a  steady  state,  which  was  considered 
to  be  the  point  at  which  surface  temperature  indications  provided  1^  tire  thermocouple 
attached  to  the  heater  sti^  stabilised.  Stabilisation  times  were  typically  between  10 
and  15  rtunutes  (this  could  be  longer  for  large  terrqrerature  changes  and/or  low  qreed 
flows).  Section  3. 1.1.1  details  a  rrrethod  by  which  better  estimates  of  stabilisation 
times  were  obtained.  Heat  was  transferred  firtan  tire  strip  via  three  medranisms-- 
cortductitm  irrto  the  perspex  plate,  radiatitm  from  die  surface  to  the  surroundings  and 
cottvective  heat  trartsfer  to  the  air  stream. 

The  geometry  of  die  plate  was  sudt  that  conduction  imo  the  per^rmi  was 
tregligible  once  steady  state  was  obtained.  This  was  achieved  widi  a  flat  plate  and 
heater  arrangement  in  whtoh  flows  on  eidier  side  of  the  plate  were  assumed 
symmetrical.  Transition  locations  were  assumed  to  be  identical  on  eidier  side  of  the 
plate.  Transverse  conduction  into  the  surrounding  per^iex  at  the  edges  of  the  heater 
str^  was  estimated  to  be  a  maximum  of  8%  (for  a  fuUy  laminar  bomidary  layer). 
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However,  results  were  only  taken  at  the  centre-line  of  the  strip,  where  estimates 
suggested  lateral  conduction  to  be  negligible.  Conservative  estimates  found 
conduction  in  the  longitudinal  direction  to  be  a  maximum  of  0.8%  of  die  convective 
heat  transfer  for  the  perspex  plate  (see  Appendix  C). 

Due  to  the  uncertainty  associated  with  the  emissivity  of  the  LC  layer  no 
radiation  correction  has  been  made.  BAUGHN  et  al.  (198S)  suggest  that  emissivity 
varies  considerably  with  the  thickness  of  the  LC  layer.  Infrared  measurements  made 
(under  the  assumption  of  a  grey  botfy)  on  a  body  similar  to  the  perspex  model,  but 
with  ink  (rather  than  matt  black  paint)  lying  over  die  crystal  layer  (rather  than  under 
it),  gave  an  emissivity  of  0.96.  This  agrees  roughly  with  typical  emissivities  given  by 
BAUGHN  et  al.  Estimates  based  on  these  emissivities  indicate  that  radiation  could 
amount  to  a  maximum  of  7%  of  the  convective  heat  transfer. 

Only  the  upper  surface  of  the  plate  was  spaifoi  with  Mack  paint,  while  die 
aluminium  was  left  exposed  on  the  lower  surface.  Therefore,  owing  to  the  relatively 
lower  radiation  on  the  lower  surfece,  a  small  temperature  gradient  will  develop 
between  the  upper  and  lower  surfaces.  However,  estimates  suggest  that  the  additional 
heat  flux  into  the  upper  surface  is  less  than  0.5%. 

Markings  were  {daced  alongside  the  heater  in  10  mm  increments  from  die 
leading  edge.  The  locations  of  liquid  crystal  isotherms  aid  gas  temperature  were 
noted  and  a  local  heat  transfer  coefBciem,  h,  was  calculated  from 
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k  = 


i*e 


oe(r/-r/) 


(2J) 


where  i,  is  the  cuirent,  w  is  die  width  of  the  heeter,  is  the  calibrated  ten^rature 
of  the  liquid  crystal,  7^  is  the  gas  temperature,  o  is  die  Stefan-Boltzman  constant,  e 
is  the  emissivity  (which  was  set  to  zero  for  the  present  tests)  and  6  is  the  resistance 
per  square  of  the  heater.  Local  O  was  corrected  for  tenqierature  variations  in  resis¬ 
tance  using, 

e  =  e„,[  1+0(7, -r„,)]  .  (2J) 

where  a  is  tenqierature  coefficient  of  resistivity  and  the  subset^  indicates  values 
at  a  reference  condition.  This  was  repeated  for  a  number  of  cuirent  settings. 
Isotherms  tqqpeaied  in  different  locadons  for  eadi  pown  setting,  thereby  allowing  h 
to  be  determined  at  a  number  of  locations  along  die  lengdi  of  the  plate.  Transition 
start  and  end  were  identified  as  described  in  Sectitm  2.2. 


2.2.1.2,2  Flat  Plate,  Constant  Pressure  Gradient 

In  many  practical  {pplications  in  which  the  steady-state  heat  transfer  technique 
enqiloying  teiiiperatuie  sensitive  LC’s  may  be  used  to  detect  transition  location  or  nup 
Stanton  number  contours,  thoe  exists  a  non-zero  }»essuie  gradi«it.  The  NLP  project 
is  paiticulariy  relevant.  The  geometry  of  the  NLF  nacelle  is  sudi  that  a  fiivourable 
pressure  gradient  exists  over  a  large  portion  of  the  beater.  Hierefore,  temperature 
sensitive  liquid  crystal  experimeiils  were  performed  under  a  favouraUe  {nessuie 
gradient  widi  the  aims  of  (1)  simulating  the  NLF  enviroiunent,  (2)  determining 
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whether  present  e]q)eriinental  techniques  are  suitable  for  detecting  transition  location 
in  a  favourable  pressure  gradient,  (3)  investigating  the  influence  of  surface  heating  on 
boundary  layer  transition  in  favourable  pressure  gradients  and  (4)  deteimining  whether 
a  Stanton  number  could  be  used  to  coiielate  the  movements  of  transition  with 
variations  in  heating  in  a  manner  similar  to  that  used  when  no  pressure  gradient  was 
present. 

Favourable  pressure  gradients  were  aj^lied  to  the  plate  using  an  adjustable  flat 
fairing  as  described  in  Section  2.1.3.  This  produced  an  approximately  unifoim 
pressure  gradient  on  the  test  surface.  An  approximately  zero  pressure  gradient  was 
achieved  wifli  die  fairing  in  the  most  upnight  position.  Details  of  the  pressure  gradient 
are  given  in  Section  3.1. 1.2. 

The  perspex  plate  and  aluminium  heater  used  fin  die  favourable  pressure 
gradiem  tests  were  the  same  as  diose  used  for  the  zero  pressure  gradient  tests. 
However,  die  means  of  gathering  data  was  slighdy  di^rem.  Previoufly,  die  iiqiut 
current  was  varied  aibitrarily  in  order  to  produce  a  series  of  constant  temperature 
contours  (refer  to  Section  2.2.1.2.1).  In  die  pressure  gradient  tests,  however,  the 
current  was  increased  until  die  LC  just  began  to  show  colour  (indicating  the  transition 
start  location).  The  current  was  then  incrementally  reduced  as  the  plate  stabilised  sudi 
that,  when  the  plate  had  reached  the  fully  stabilised  cmditirai,  the  surface  was  being 
heated  just  enough  to  provide  an  indication  of  transition  start.  Similatly,  the  transition 
end  location  and  its  associated  heating  rate  were  established  by  varying  die  current 
sudi  that,  when  die  pltue  was  stabilised,  the  contours  of  coratant  temperature  moving 


downsoeam  and  those  moving  upstream  just  touched.  The  process  was  repeated  for 
each  of  the  three  LC’s  at  eadi  fairing  setting.  This  method  gave  three  transition  start 
and  end  locations  and  the  heating  rate  associated  widi  each  for  a  range  of  pressure 
gradients.  Previously,  these  data  had  to  be  interpolated  from  plots  of  h  versus  x  such 
as  shown  in  Fig.  2.3a.  and  were,  therefore,  less  accurate. 

Since  the  plate  with  the  aluminium  heating  element  did  not  contain  pressure 
ttqjpings,  it  was  not  possible  to  diiecdy  evaluate  the  local  pressure  gradient  and  flow 
speed  when  the  aluminium  heater  was  in  the  tunnel.  Therefore,  the  plate  with  static 
pressure  tilings  was  used  to  directly  measure  the  pressure  gradirats  and  flow  speeds. 
These  results  were  noted  and  related  to  the  heated  plate  through  the  oncoming  flow 
speed  and  pressure  fairing  setting  as  described  in  the  next  paragr^h. 

With  the  pressure  plate  in  the  tunnel,  the  pressure  fairing  was  run  through  a 
series  of  settings  labelled  settings  0-10.  Setting  zero  corresponded  to  the  fairing  being 
drawn  fully  up  (zero  pressure  gradient).  At  each  successive  setting,  the  downstream 
support  was  lowered  5  mm.  The  upstream  support  was  adjusted  so  as  to  make  the 
pressure  gradient  as  uniform  as  possible.  The  positicm  of  each  support,  the  oncoming 
speed,  total  pressure  and  local  static  pressures  were  then  recorded.  This  procedure 
was  rqteated  for  each  tunnel  speed  setting  of  interest.  The  positions  of  the  supports 
which  provided  the  most  uniform  pressure  gradkot  for  a  given  fairing  setting  were 
found  to  be  independent  of  turmel  ^leed  setting.  At  each  fairing  setting,  local 
velocities  on  the  surface  of  die  plate  were  related  to  oncoming  speed  by  fitting  a  first 
order  least-squares  curve  fit  to  a  {riot  of  distance  along  die  plate,  x,  versus  the  ratio 
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of  local  to  oncoming  velocities.  Hence,  when  the  heated  plate  was  installed  in  the 
tunnel,  die  oncoming  speed  was  measured  directly  and  local  speed  was  found  from 

u,  =  uiA*Bx)  ,  (2.4) 

where  u,  is  the  local  velocity,  u  is  the  velocity  at  the  leading  edge  of  the  plate 
(oncoming  velocity  corrected  according  to  Eqn  2.1),  and  A  and  B  are  constants. 
Values  of  A  and  B  for  all  conditions  of  interest  to  these  tests,  as  well  as  the 
corresponding  positions  of  both  supports,  are  provided  in  Table  2.1. 


Table  2.1  Pressure  Fairing  Schedule 


Fairing 

Setting 

Support  1 
(mm) 

Suf^rt  2 
(mm) 

Oncoming 

Speed 

(m/s) 

A 

B 

0 

0 

0 

:»).s 

1.00 

0.00 

1 

2 

4 

20.7 

1.00 

0.02 

2 

4 

9 

20.5 

0.99 

0.04 

0 

0 

0 

26.4 

1.00 

0.00 

1 

2 

4 

26.3 

1.00 

0.01 

2 

4 

9 

26.1 

0.99 

0.03 

3 

5 

14 

25.8 

0.98 

0.06 

0 

0 

0 

32.4 

1.00 

0.01 

1 

2 

4 

34.2 

0.94 

0.02 

2 

4 

9 

31.9 

0.99 

0.04 

3 

5 

14 

31.6 

0.98 

0.06 

4 

9 

19 

31.2 

0.98 

0.09 

Note  that  the  oncoming  speed  tended  to  decrease  as  the  pressure  fairing  was 
lowered,  contracting  the  test  section  and  diereby  increasing  the  drag  on  die  motor 
(which  was  held  at  a  constant  power  setting).  The  constants  A  and  B  are  dependent 
on  failing  setting,  but  show  litde  dependence  on  oncoming  flow  speed.  Neveitheless, 
A  and  B  were  only  employed  when  the  oncoming  flow  over  the  heated  plate  was 
within  2  tn/s  of  that  for  which  the  constants  were  determined.  This  led  to  estimated 
errors  in  local  velocity  of  <  1%. 

2.2  J,  Carbon  Heater 

2.2,2.1  Experimental  Apparatus 

Oxford  University  is  to  perform  in-flight  heat  transfer  experiments  on  the 
surface  of  an  NLF  nacelle.  Experimental  heat  transfer  distributions  will  be  used  to 
determine  the  location  of  transition  from  laminar  to  turbulent  flow  for  a  range  of  test 
conditions.  A  heater  is  to  be  incorporated  into  the  skin  of  the  nacelle  and  isolated 
firom  the  surrounding  carbtm  structure.  The  heater  wfll  consist  of  a  str4>  of 
unidirectional  carbon  fibres  with  uniform  resistance  along  its  loogth.  A  uniform 
surface  heat  flux  will  be  created  by  passing  a  current  thiou^  the  fibres.  The  present 
experiments  are,  to  die  audior’s  knowledge,  the  first  to  attempt  to  locate  transition  and 
make  quantitative  heat  transfer  measurements  employing  part  of  the  aircraft  strocture 
(the  carbon  fibre)  as  a  heating  element. 

A  flat  model  of  the  carbon  heater  which  is  to  be  incorporated  imo  the  NLF 
nacelle  was  manufactured  to  the  specifications  of  Oxford  University  by  Hurel-Dubois 
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(UK)  ltd.,  of  Burnley,  England.  Incon^piessible  steady  state  convective  heat  transfer 
tests  similar  to  those  performed  cm  the  aluminium  foil  were  carried  out  on  die  carbon 
sanqile.  These  experiments  were  intended  to  test  the  thermal  and  mechanical 
properties  of  the  carbon  satiqile,  and  to  ensure  that  the  methods  of  detecting  transition 
which  were  developed  on  the  aluminium  plate  remained  valid  for  the  carbon  sample. 

The  upper  surface  of  die  carbon  structure  is  shown  mounted  in  die  perspex 
plate  in  Figure  2.9.  The  carbon  model  measured  0.203  m  x  0.518  m,  with  a  0.100  m 
X  0.565  m  heating  element  centred  in  the  model.  The  model  was  mounted  in  a 
perspex  plate  identical  to  die  one  described  in  Section  2.1.2  with  the  heating  element 
down  the  centre-line  of  the  plate.  After  smoothing  with  plastic  padding,  the  step  from 
the  perspex  to  the  carbon  plate  was  <  50  pm  at  a  position  65  mm  downstream  of  the 
leading  edge.  The  unheated  staiting  loigdi  was  115  mm.  There  was  no  detectable 
step  from  die  heater  to  the  surrounding  carbon-fibre  stiucture. 

The  heating  elemem  consisted  of  a  single  layer  of  high  tensile  strength  (HTS) 
unidirectional  PAN  fibres  widi  a  nominal  dikkness  of  0.125  mm.  The  unidirectional 
fibres  were  insulated  on  the  underside  by  two  liters  of  0.3  mm  woven  glass  cloth 
(weave  style  7781)  with  a  diermal  conductivity,  k,  of  0.03  W/(m  K)  (refer  to 
TOULOUKIAN  et  al.  (1S170)).  Both  the  beater  and  glass  cloth  were  set  into  a  woven 
carbon  cloth.  The  entire  assembly  was  inqnegnated  with  epoxy  resin  which  is  also 
an  insulator,  k  «>  0.23.  The  nominal  diickness  of  die  con^iosite  was  ~  3  mm,  and  the 
fibre  volume  firaction  was  57  %  (reference  CLAYTON  (1991)). 
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The  unidiiectional  caibon  fibres  which  composed  the  heating  strip  had  a 
unifonn  resistance  in  the  longitudinal  direction  of  0.22  Q/square.  A  qualitative 
measure  of  die  uiuformity  of  the  heater’s  resistance  was  made  by  passing  a  current 
through  the  heater  under  zero  flow  conditions  and  checking  for  uniformity  in  the 
colour  play  of  the  LC’s  on  its  surface.  Resistance  was  found  to  be  extremely  uniform. 
The  carbon  heater’s  temperature  coefficient  of  resistance,  a,  was  measured  to  be  -  0.0 
by  placing  a  sample  of  the  carbon  fibre  in  a  bath  of  de-ionised  water  and  measuring 
its  resistance  over  a  temperature  range  of  22-70°C. 

Three  temperature  sensitive  liquid  crystals  were  applied  to  the  surface  of  the 
caibon  structure  in  the  manner  which  was  detailed  in  Section  2.2.1. 1.  Figure  2.10 
shows  the  complete  packaging  arrangement  of  the  heater  and  LC’s.  As  the  surface 
of  the  carbon  composite  is  already  black,  painting  was  unnecessary.  A  type  K 
thermocouple  of  125  pm  thickness  was  mounted  on  the  surface  of  the  heater  for  use 
in  calibrating  the  LC’s. 

RHODES  and  MOSES  (1991)  have  noted  that,  due  to  its  anisotropic  nature, 
both  thermal  and  electrical  resistance  in  unidirectional  carbon  fibres  are  higher  in  the 
transverse  direction  than  in  the  longitudinal  direction  by  a  factor  of  ^iproximately  10. 
Therefore,  poor  electrical  contact  was  found  to  produce  non-uniform  heating  on  the 
test  surface,  in  addition  to  contact  heating.  In  order  to  allow  for  electrical  contact 
without  disturbing  the  finish  on  the  upper  surface,  the  unidirectional  fibres  were  fed 
through  a  slit  in  the  insulating  glass  cloth  and  in  the  woven  carbon  cloth  and  extended 
out  the  bottom  of  the  stmctuie.  Fig.  2.11  shows  the  underside  of  the  carbon 
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con^site  mounted  in  the  peispex  plate.  The  best  electrical  connection  was  made 
aiiter  removing  the  resin  from  the  ends  of  the  fibres.  These  were  then  clanqied  onto 
a  metal  frame  which  was  glued  to  the  underside  of  the  carbon-fibre  stmcture.  The 
free  carbon  fibres  were  clatrqxd  wifb  a  str^  of  mbber  which  deformed,  allowing  all 
fibres  to  come  into  contact  with  the  metal  base.  Fig.  2.12  shows  a  cross-section  of 
the  electrical  contact.  A  photo  of  the  contact  is  provided  in  Figure  2.13a. 

It  is  desirable  to  minimise  the  length  of  fibre  which  extends  from  the  heater 
to  the  electrical  contact.  Upon  leaving  the  connection,  the  fibres  are  generating  heat. 
This  could  lead  to  overheating  in  the  nacelle  if  the  heat  is  not  dissq>ated.  Therefore, 
the  amount  of  heat  generated  along  the  length  of  the  carbon  fibre  leads  was  reduced 
by  soaking  the  fibres  in  electrically  conducting  paint,  see  Figure  2.13b.  This 
effectively  increased  the  length  of  the  contacts  down  to  die  point  where  the  fibres 
remain  coated  with  epoxy. 

A  silver  loaded  conducting  paint  based  on  acetone  and  xylene  was  chosen  to 
reduce  the  resistance  of  the  free  fibres.  The  paim  was  supplied  by  RS  Conqionents 
of  London,  England.  The  fibres  were  soaked  in  the  conducting  paint  which  was 
diluted  with  acetone  in  order  to  allow  the  paint  to  better  penetrate  between  the  fibres. 
The  acetone  then  evrqxirated,  leaving  the  conducting  paint. 

The  effective  reductitm  in  overall  heater  length  was  i^arent  from  the 
reduction  in  the  carbon  heater’s  total  resistance.  The  silver  loaded  paint  was  applied 
to  -  4S  mm  of  die  fibre  leads  and  reduced  the  heater’s  total  resistance  by  0.10  fl. 
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Thus,  given  a  heater  resistance  of  0.22  fVsquaie,  the  heater’s  length  was  effectively 
reduced  by  the  length  of  the  carbon  diat  was  treated  with  the  conducting  paint. 

2J!,2.2  Experimental  Program 

Current  was  passed  through  the  heating  str^,  creating  a  urufotm  surface  heat 
flux.  As  with  the  perspex  plate,  heat  was  transferred  from  die  heated  surface  via 
radiation  firom  the  surface  to  the  surroundings,  conduction  to  the  surrounding  carbon 
structure  and  convection  to  the  air  stream.  Radiation  was  ignored  due  to  the  uncertain 
emissivity  of  the  LC  layer  (refer  to  Section  2.2.1. 2.1). 

An  effort  was  made  to  minimise  the  thermal  conductivity  of  the  carbon 
stmcture  through  careful  selection  of  die  carbon  fibres  from  which  it  was  constmcted. 
Nevertheless,  thermal  conduction  was  found  to  be  much  mote  significant  in  die  carbon 
model  than  in  the  perspex  model. 

A  number  of  studies  have  investigated  the  tbential  conductivity  of  unidirec¬ 
tional  carbon  fibres  in  both  the  longitudinal  and  transverse  directions.  The  thermal 
conductivity  of  some  woven  fabrics,  such  as  Kevlar,  has  also  been  studied  (TOULOU- 
KI.\N  et  al.  (1970)).  However,  to  die  author’s  knowledge  there  are  no  studies 
addressing  the  conductivity  of  woven  carbon  fabrics.  Therefore,  a  model  was 
constmcted  which  assumed  that  conductioo  in  the  woven  carbon  fibres  occurred  solely 
in  the  longimdinal  direction  of  die  fibres.  As  the  fibres  in  the  woven  cloth  run  at  45 
degree  angles  to  the  heater’s  unidirectional  fibres,  this  effectively  reduced  die 
longitudinal  thermal  conductivity,  k,  by  l/>/2  in  the  direction  perpendicular  to  die 
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heater.*  The  longitudinal  thennal  conductivity  of  the  fibres  was  estimated  to  be  9.5 
W/(m-K)  at  BOO^K  based  upon  die  estimates  of  CLAYTON  (1990)  as  well  as  die 
woik  of  AGARWAL  and  BROUTMAN  (1980)  and  RHODES  and  MOSES  (1991). 

Lateral  thermal  conduction  from  the  heater  to  the  surrounding  carbon  structure 
was  found  to  be  significant  in  relation  to  convective  heat  transfer.  However, 
measurements  were  taken  only  on  the  centre-line  of  the  heater,  where  errors  due  to 
conductive  losses  were  estimated  to  be  2-8%  of  the  convective  heat  transfer  rate.  This 
estimate  assumed  a  linear  tenqwrature  gradient  from  the  centre-line  to  the  edges  of  the 
heater.  In  reality,  the  temperature  gradient  is  nearly  zero  at  the  centre-line  of  the 
heater.  Hence,  lateral  conductive  losses  wiU  be  much  less  near  the  centre-line. 
Conductive  losses  will  be  further  reduced  on  the  NLF  nacelle  where  the  thermal 
conductivity  of  the  carbon  stracture  will  be  reduced  by  -  20%  due  to  the  reduced 
thermal  conductivity  of  unidirectional  fibres  at  lower  ten^)eratures.  Fig.  2.14  shows 
the  variation  of  thermal  conductivity  witii  tenqterature  for  unidirectional  fibres.  This 
plot  was  constmcted  from  die  results  of  RHODES  AND  MOSES  (1991). 

The  underside  of  the  carbon  plate  was  insulated  with  25  mm  of  styrofoam. 
Conduction  to  the  underside  of  the  plate  was,  therefore,  considered  negligible.  It  was 


*  Note  that  longitudinal  tiiennal  conductivity  of  the  carbtm  plate  is  largest  for 
the  case  in  which  the  fibres  witiiin  the  layered  cloth  run  at  45°  angles  to  the 
unidirectional  fibres.  As  the  carbon  cloth  is  rotated  in  relation  to  the  unidirectional 
fibres,  the  effective  k  of  the  cloth  is  continuously  reduced.  When  die  cloth  has  been 
rotated  a  full  90°,  so  as  to  {dace  successive  layers  of  die  woven  carbon  fibre  paraUel 
and  perpendicular  to  the  unidirectional  fibres,  the  L  of  the  carbon  structure  in  die 
longitudinal  direction  has  been  effectively  reduced  to  1/2  the  value  of  the  Itmgitudinal 
thermal  conductivity  of  die  carbon  fibres. 
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also  found  to  be  in^Mitant  to  insulate  the  contacts  and  the  area  surrounding  diem.  If 
left  uninsulated,  heat  was  rtpidly  convected  away  from  this  area,  paiticulariy  as  the 
£ree-stieam  air  flowed  through  the  exposed  cadxm  fibre  leads.  This  resulted  in 
conduction  from  other  sections  of  the  plate,  producing  spurious  heat  transfer  results. 

A  worst  case  analysis  of  longitudinal  conduction  is  presented  in  ^ipendix  C. 
The  regions  where  errors  due  to  conduction  ate  expected  to  be  the  greatest  are  diose 
in  which  is  at  a  maximum  or  minimum.  These  regions  correspond  roug^y 

to  the  Stan  and  end  of  transition.  Results  indicate  that  the  calculated  values  for 
Stanton  number  at  the  stan  of  transition  may  be  as  much  as  13%  high.  Again,  diis 
may  be  reduced  to  ~  11%  for  the  NLF  nacelle  due  to  the  drop  in  die  theimal 
conductivity  of  the  carbon  fibre  widt  decreasing  temperature  (refer  to  Fig.  2.14). 
Conversely,  at  transition  end,  calculated  values  of  Stanton  number  may  be  as  much 
as  13%  percent  too  low  on  the  carbon  model.  Theie  will  be  litde  error  in  areas  where 
the  tenqieratuie  gradimt  is  small  (sudi  as  in  the  turbulent  region).  Conduction  does 
not  affect  the  location  of  transition  as  indicated  by  die  LC’s.  However,  as  noted  in 
.^ipendix  C,  the  calculated  location  of  adiabatic  Reynolds  number  may  be  slighdy 
affected. 

Current  was  supplied  to  the  heater  and  the  carbon  model  was  allowed  to  reach 
thermal  stability  (typically  taking  on  the  order  of  4  minutes).  The  distribution  of  heat 
transfer  coefficients  along  the  length  of  the  plate  was  found  using  die  method 
described  in  Section  2.2.I.2.I. 
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23  Shear  Sensitivt  Liquid  Crystal  Experiments 

Experiments  to  measure  surface  shear  stress  were  perfoimed  using  a  shear 
sensitive,  tenqterature  insensitive,  unencqMulated  diolesteric  liquid  ciystal.  Previous 
uses  of  shear  sensitive  crystals  for  aerodynamic  studies  (e.g.  HOLMES  et  al.  (1986) 
and  REDA  (1991))  have  enq>loyed  the  reversible  colour  play  technique  discussed  in 
Section  1.3.  This  study  represents  the  first  experimental  study  implementing  the 
irreversible  texture  change  method,  first  introduced  by  BONNETT  (1989),  to 
measurements  in  transitional  boundary  layers.  Using  this  technique,  the  LC  is  applied 
to  the  test  surface  and  is  drivm  into  the  coloudess  focal-conic  texture.  If  this  is  dren 
sheared,  the  movement  of  the  film  aligns  the  cholesteric  helices  and  the  layer  displays 
the  colourful  Grandjean  texture  (see  BONNETT  et  al.  (1989)).  This  process  of 
changing  texture  is  irreversible  unless  the  LC  is  heated  to  drive  it  back  into  the  focal 
conic  texture.  However,  unlike  the  convoitional  reversible  technique,  there  is  no 
dependence  on  illumination  or  viewing  angle.  Hierefore,  fix  LC  may  be  calibrated 
to  give  a  quantitative  measure  of  shear  stress. 

Z3.1  Experimental  .^^aratus 

A  flat  plate,  similar  to  that  used  for  beat  transfer  measurements,  was  utilised 
for  the  shear  stress  measurements.  The  back  surface  of  the  plate  was  qrrayed  vrith  a 
matt  black  paint  to  clarify  the  colour  change  of  the  LC  li^er.  The  iteat  diolesteric 
liquid  crystal,  designated  HSll,  was  siq)plied  by  BDH  Ltd.  of  Poole,  England. 

The  cholesterics  can  be  dqiosited  on  the  surface  by  spraying,  silk-screening  or 
painting.  Using  the  latter  technique  it  was  difficult  to  obtain  a  film  of  uniform 
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thickness.  In  order  to  spray  the  coating  <Mito  the  surface  the  liquid  crystal  is  first 
dissolved  in  a  suitable  solvent,  typically  trichloroethane.  After  spn^ing,  the  solvent 
evqMralBS  leaving  a  thin  film  of  cholesteric  LC  on  the  surface.  This  technique  of 
q>plication  has  been  used  successfully  by  REDA  (1991)  and  HOLMES  et  al.  (1986) 
in  aerodynamic  flight  tests  in  which  a  thin  coating  (-  10pm)  of  LC  was  desirable, 
such  that  the  film  changed  colour  but  did  not  run  under  aerodynamic  loading.  The 
texture  change  technique,  however,  requires  that  the  LC  layer  be  thick  enough  to  allow 
the  crystal  layer  to  flow.  This  thickness  was  difficult  to  achieve  when  spn^ing. 
Furthermore,  when  sprayed  on,  each  coating  could  be  used  only  once  and  it  was 
difficult  to  ensure  that  all  of  the  solvent  had  evtporated. 

Silk-screening  was  found  to  be  the  best  technique  for  the  present  experiments. 
This  produced  a  film  of  uniform  thickness  which  could  be  re-used  by  re-screening  the 
same  mixture  after  it  had  flowed  during  an  experiment.  The  results  presented  in 
Section  3.2  were  obtained  widi  a  silk-screened  layer  q)proximately  40  pm  thick. 

It  was  found  to  be  iitqxmam  to  avoid  dust  contamination  of  the  surface.  Dust 
particles  on  the  LC  film  sometimes  moved  under  the  action  of  aerodyiumic  drag 
resulting  in  the  formation  of  streaks  of  Grandjean  texture  in  the  film.  Any  inqterfec- 
tions  or  bubbles  in  the  LC  film  were  manually  tentoved,  and  die  LC  was  re-screened. 

The  cholesteric  LC  was  dqwsiied  along  a  str^  60  mm  wide  and  700  mm  Irmg 
down  the  centre  of  the  plate.  The  film  started  60  mm  fiom  the  leading  edge  of  the 
plate.  This  was  done  to  prevem  the  LC  layer  from  disturbing  the  boundary  Ir^er  at 
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the  firont  of  the  plate  tK^tere  the  boundaiy  layer  was  thin  and  thereby  inducing  early 
transidon.  Calculations  based  on  die  review  of  YOUNG  and  PATTERSON  (1981) 
indicate  that  a  two  diniensitMul  excrescence  of  bei^  less  than  70  pm  at  60  mm  from 
the  leading  edge  will  not  distuib  the  boundary  layer  for  any  of  the  condititnis  of  die 
present  tests. 

J,  ExperimeiUal  Program 

After  silk  screening,  most  of  the  LC  layer  was  in  the  Grandjean  texture  and  the 
layer  was  a  dull  red  colour.  The  film  was  driven  back  to  the  colourless  focal  conic 
state  by  heating  it  into  the  isotropic  phase  (approximately  50°C  for  the  71511  mixture 
used)  with  a  hot  air  gun  and  allowing  it  to  cool  back  to  the  cholesteric  phase. 

The  tunnel  was  staited  when  the  plate  had  letumed  to  ambient  temperature, 
applying  steady  shear  stress  to  the  measurement  surface.  After  die  shear  has  been 
applied,  it  takes  some  time  fOT  die  film  to  show  a  duU  red  colour.  As  the  cholesteric 
helices  align  with  increased  displacement  of  die  film  surface  (i.e.  as  the  mesophase 
progresses  from  the  focal-conic  to  the  Grandjean  texture)  the  viscosity  of  the  layer 
rapidly  decreases.  Eventually  the  liquid  crystal  film  flows  and  displays  a  de^  red 
colour. 


Markings  were  placed  alongside  die  LC  film  at  10  mm  increments  from  the 
leading  edge.  The  location  at  which  the  LC  began  to  flow  aiKl  the  length  of  time  for 
which  the  shear  stress  had  been  aq^died  were  noted.  The  location  at  which  the  LC 
film  first  begins  to  show  colour  (region  of  highest  shear  stress)  provides  a  qualitative 
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indication  of  transition  end.  With  time,  the  location  at  which  the  ciystal  begins  to 
flow  moves  forward  throu^  the  transitional  region,  and  downstream  through  the 
turbulent  region.  The  time  titat  is  required  to  induce  flowing  in  the  LC  film  is  used 
to  quantify  dte  surface  shear  stress  as  discussed  in  Section  3.2. 
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Chapter  3:  Results  and  Discussion 


3.1  Temperature  Sensitive  Liquid  Crystals 

3.1.1  Aluminium  Heater 

3.1.1.1  Flat  Plate,  Zero  Pressure  Gradient 

Basic  results  from  the  heat  transfer  experiments  are  presented  for  a  range  of 
flow  speeds  in  Fig  3.1.  Results  are  shown  for  all  three  liquid  crystals  at  each 
conditicMi  and,  for  conqrarison,  laminar  and  turbulent  predictions  from  the  well 
established  boundary  layer  code  STANS  (CRAWFORD  and  KAYS  (1976)).  In  the 
regions  of  fully  laminar  and  fully  turbulent  flow,  the  levels  of  Stanton  number 
indicated  by  each  of  the  crystals  are  in  good  agteetnem.  In  general,  the  levels  of  heat 
transfer  ate  lower  than  would  be  »pected  from  standard  conelations  and  predictions 
of  STANS,  particularly  in  the  tegitm  of  fully  turbulent  flow.  This  is  also  apparent  in 
the  results  of  KIM  et  al.  (1989)  in  fiilly  turbulent  regions. 

Experimental  uncertainty  alone  does  not  account  for  the  low  levels  of  heat 
transfer  in  the  turbulent  region.  Uncertainties  in  measured  values  of  heat  transier 
coefficient  are  estimated  to  be  maximum  of  ±  8.7%  in  die  turbulent  region  for  these 
experiments.  There  are  at  least  two  other  possible  explanations  for  unexpectedly  low 
measurements  of  the  coefficient  of  heat  transfer.  Rrst,  as  mendtHied  in  Section  2.1.1, 
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flow  visualisation  indicated  that  there  may  be  smne  sectmdaiy  flows  generated  at  the 
tunnel  contraction  resulting  in  a  non-two-dimensional  boundary  layer  on  the  plate. 
Therefore,  a  section  of  the  heater  may  have  experienced  a  relatively  low  mass  flow 
rate  over  its  sur^e.  The  second  possibility  is  that  an  area  of  relatively  high 
resistance  existed  within  the  heater.  Either  possibility  would  cause  a  local  ’hot  spot’ 
to  develop.  As  a  result,  the  local  coefficient  of  heat  transfer  would  be  unexpectedly 
low.  Such  a  ’hot  spot’  was  apparent  on  the  heating  strip  at  -  3S0-SS0  mm  from  the 
leading  edge.  Nevertheless,  the  results  suggest  that  the  arrangement  was  suitable  for 
the  present  purposes. 

The  results  clearly  show  the  influence  of  heating  on  transition  location.  The 
lowest  temperature  LC  consistently  shows  a  later  start  and  end  of  transition.  This  is 
because  the  lower  teirqterature  IX  displays  colour  at  a  lower  heating  rate  than  the 
higher  temperature  LX’s.  As  the  heating  rate  is  increased,  both  start  and  end  of 
transition  are  moved  forward. 

The  shape  of  the  heat  transfer  rise  through  die  transitional  region  and  the 
length  of  transition  are  distorted  when  presented  as  in  Fig.  3.1.  For  a  given  LC,  the 
heating  requited  to  produce  colour  change  at  the  start  of  transition  is  less  than  that 
required  to  produce  change  at  the  end  of  transition.  Thus,  when  the  Stanton  number 
is  being  measured  at  the  end  of  transition,  the  transition  start  location  will  have  moved 
further  upstream  because  of  the  forward  movement  of  transition  with  increased 
heating.  It  is  possible,  therefore,  to  see  a  maximum  Stanton  number  upstream  of  die 
minimum  in  these  plots. 


3.2 


The  2q^>arent  length  of  transition  can  also  appear  distorted.  For  example,  when 
coefficient  of  heat  transfer  is  plotted  so  as  to  indicate  the  heat  transfer  distribution 
along  the  plate  as  in  Fig.  3.2,  the  length  of  transition  q^iears  misleadingly  short.  At 
any  given  heating  rate,  the  start  and  end  of  transititHi  would  be  further  ^art  and  the 
rise  through  transition  mote  gradual  than  is  indicated  in  these  figures. 

As  previously  mentioned  (Section  1. 2.2.4),  it  is  well  established  that  a  surface 
heat  flux  destabilises  a  laminar  boundary  layer  leading  to  early  transition.  The 
magnitude  of  this  effect  for  a  given  flow  and  heat  flux  is  of  importance  if  die 
technique  is  to  be  used  in  an  attempt  to  determine  where  transition  would  occur  on 
an  unheated  sutface  (as  is  its  intended  use  for  the  NLF  program).  By  using  three 
liquid  crystals  on  the  heater  simultaneously,  this  effect  was  examined  by  observing  the 
change  in  transition  location  with  different  surface  heat  fluxes. 


Dimensional  analysis  (see  ^^ndix  A)  for  fliis  inconqiressible  flow  indicates 
that  the  Reynolds  numbers  at  start  and  end  of  transition  should  be  a  function  of 
Stanton  number  in  the  form 


4 

p«c,r. 


(3.1) 


where  q  is  the  convective  heat  transfer  rate,  p,  u  and  are  the  density,  speed  and 
temperature  of  the  fiee-stream  fluid,  respectively,  and  C,  is  die  specific  heat  of  the 
fluid.  From  the  experimental  results,  the  Reynolds  numbers  at  the  start  and  end  of 
transition  for  each  LC  at  each  flow  speed  were  determined.  This  gave  transition 
R  .ynolds  numbers  at  three  heating  rates  at  each  qieed.  An  equation  for  critical 
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Reynolds  number.  Re^  as  a  function  of  the  heating  parameter  of  the  form 


Re 


crlt 


-B 

=  At 


i 

P.c,r, 


(3^) 


was  least-squares  fitted  to  each  set  of  three  points,  and  an  extnqjoladon  to  zero  heat 
transfer  rate  was  used  to  determine  an  adiabatic  transition  Reynolds  number,  Re^,^, 
for  each  speed.  The  results  for  both  start  and  end  of  transidcm  at  all  speeds  ate 
plotted  in  Pig.  3.3.  These  suggest  that  the  Stanton  number  parameter  used  correlates 
the  data  quite  well.  The  reduction  in  transition  Reynolds  number  is  shown  to  be  an 
exponential  function  of  Stanton  number.  The  maximum  heating  effect  is  rou^y  a 
halving  of  the  adiabatic  critical  Reynolds  number  at  Stanton  numbers  greater  than  - 
0.3  X  10^.  At  a  given  speed,  transition  end  qrpears  to  be  less  effected  by  the  heating 
than  is  transition  start.  However,  the  difference  is  slight  and  the  effect  on  both 
transition  end  and  stan  may  reasonably  be  related  by  a  single  curve.  A  least-squares 
curve  fft  to  the  data  has  the  form 

Du  -10.4x10* — 1— 

- =  1  -0.50(l-e  _  (33) 

^erU^ 


The  results  shown  in  Fig.  3.3  agree  qualitatively  with  the  experimental  results 
of  LAUCHLE  and  GURNEY  (1984).  Their  plot  of  transition  Reynolds  number  versus 
heating  (Fig.  1.16)  is,  qualitatively,  a  minor  image  of  the  results  obtained  in  this 
experiment.  This  is  to  be  expected  since  their  experiments  were  carried  out  in  water 
(refer  to  Section  1 .2.2). 


The  time  required  for  the  per^x  plate  to  teach  thermal  stability  was  estimated 
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by  setting  the  power  which  was  being  injnit  to  the  heater  and  observing  how  die 
ten^rature  on  the  surface  of  the  plate  varied  widi  time.  The  locations  of  LC 
contours  were  noted  every  60  seconds.  For  the  first  power  setting,  the  plate  was 
heated  from  a  cooled  condition.  Subsequent  readings  were  taken  by  changing  die 
power  setting  on  the  already  heated  plate.  Calculated  heat  transfer  distributions  for 
a  transitional  boundary  layer  at  times  ranging  from  3-20  minutes  ate  shown  in  Fig. 
3.4.  The  results  suggest  diat  -  IS  minutes  are  required  for  the  perspex  plate  to  teach 
full  thermal  equilibrium.  Five  minutes  are  adequate  for  obtaining  both  qualitative  and 
quantitative  results  (to  widiin  -  2%)  in  regions  of  high  tenqieratute  gradients,  where 
a  relatively  large  change  in  die  plate’s  surface  temperature  results  in  small  changes  to 
the  stream-wise  location  of  the  LC  contour. 

In  the  turbulent  boundary  layer,  where  dTIdx  is  modest,  a  relatively  small 
change  in  surface  temperature  moves  the  LC  contours  a  significant  distance. 
However,  dh/dx  is  also  small  in  the  turbulent  regions,  such  that  the  resulting  errors  in 
the  calculated  level  of  heat  transfer  ate  small  (within  -  8%)  after  10  minutes. 

3.1.1.2  Flat  Plate,  Constant  Pressure  Gradient 

Favourable  pressure  gradients  were  iqiplied  to  the  plate  using  the  adjustable 
fairing  described  in  Section  2.1.3.  This  produced  an  qiproximately  uniform  pressure 
gradient  on  the  test  surface.  An  approximately  zero  pressure  gradient  was  obtained 
with  the  fairing  in  its  most  uptight  position.  A  typical  series  of  pressure  distributions 
measured  on  the  surface  for  various  fairing  settings  is  presented  in  Fig.  33  where 
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(3.4) 


Here,  Pg  is  the  total  pressure,  is  the  static  pressure  at  a  given  distance  along  the 
plate  and  is  a  reference  static  pressure.  The  static  pressure  at  the  first  taping  was 
used  for  P^  for  these  tests. 


The  results  of  tests  inte:''led  to  assess  the  repeatability  of  the  pressure  gradients 
are  shown  in  Figs.  3.6  and  3.7  for  (nominal)  flow  speeds  of  35.0  and  28.7  m/s 
respectively.  Each  figure  shows  a  comparison  of  two  tests  performed  on  different 
days,  but  at  the  same  nominal  speed.  TTie  results  show  excellent  repeatability. 


Dimensional  analysis  (Aj^ndix  A)  suggests  that  the  appropriate  pressure 
gradient  parameter  for  the  present  experiments  is  the  acceleration  parameter,  K,  where 


t  / 


dUf 

dx 


(3.5) 


Here,  V  is  the  kinematic  viscosity  and  u,  is  the  local  velocity.  This  acceleration 
parameter  has  been  used  previously  by  a  number  of  researchers  including  BLAIR 
(1982),  JONES  and  LAUNDER  (1972),  and  GAUGLER  (1984). 


The  conditions  for  the  experiments  are  summarised  in  Table  3.1.  The 
oncoming  speed  is  the  flow  speed  at  the  start  of  the  test  plate.  At  higher  values  of 
K  and  at  lower  flow  speeds  transition  measurements  could  not  be  made  on  the  test 
surface.  While  these  acceleration  parameters  are  small,  they  are  of  the  same  order  as 
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Table  3 . 1  Favourable  pressure  gradient  test  conditions . 


Flow  Condi¬ 
tion 

Oncoiaing 
Speed  (m/s) 

K  X  10* 

Air  Teo^- 
erature,  T 
(“C) 

la 

22.3 

0.2 

17 

lb 

22.3 

1.1 

17 

Ic 

21.9 

2.6 

17 

2a 

28.7 

0.0 

13 

2b 

28.4 

0.6 

13 

2c 

27.8 

1.9 

13 

2d 

27.3 

3.5 

13 

3a 

35.0 

0.2 

16 

3b 

35.0 

0.8 

16 

3c 

34.1 

2.7 

14 

3d 

33.3 

4.1 

11 

those  expeaed  on  the  NUF  nacelle. 

At  a  given  condition,  there  was  a  variation  in  K  along  the  length  of  the  plate 
due  to  U,  increasing  with  distance  from  the  leading  edge.  However,  the  variation  was 
small  (a  maximum  of  3.7%  at  condition  3d),  and  assuming  AT  to  be  constant  at  each 
condition  provided  a  consistent  set  of  results.  The  pressure  fairing  could  have  been 
designed  to  create  a  constant  K  rather  than  a  constant  dPtdx,  but  this  would  have 
increased  the  conqrlexity  of  the  testing. 

At  each  test  condition,  the  Reynolds  numbers  at  the  start  and  end  of  transition 
(based  on  the  local  flow  speed)  were  determined  for  three  different  heating  rates. 
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TTiese  weie  used  to  extrapolate  to  zero  heat  flux  in  order  to  obtain  the  adiabatic 
transition  Reynolds  numbers,  Re^^.  The  method  is  described  in  Section  3.1. 1.1. 

Fig.  3.8  shows  values  of  RCcriiai  normalised  with  the  value  of  at  zero 
pressure  gradient  (1?  =  0)  and  plotted  as  a  function  of  K.  The  results  for  both  start 
and  end  of  transition  fall  on  a  single  curve  and  show  a  strong  dependence  of  transition 
location  on  pressure  gradient  in  the  present  arrangement.  An  exponential  curve  of  the 
form, 

- -  =  1  +  1.07(e®”>“°'**  -1)  ,  (3.6) 

fits  the  data  weU.  The  fact  that  the  data  fall  neatly  onto  a  single  curve  provides 
conHdence  in  the  method  of  obtaining  adiabatic  transition  locations. 

Figs.  3.9-3.13  show  Re^  normalised  with  Re„i,^  and  plotted  as  a  function  of 
Stanton  number  for  each  flow  condition.  Plots  are  provided  for  both  transition  start 
and  transition  end  at  each  condition  (except  condition  3,  where  transition  end  did  not 
occur  on  the  measurement  surface).  Stuface  heating  is  shown  to  have  a  relatively 
greater  destabilising  influence  upon  the  normalised  transition  Reynolds  number  as  the 
pressure  gradient  becomes  increasingly  favourable. 

In  Section  3. 1.1.1,  it  was  shown  that  Re^Re^,^  could  be  correlated  with  a 
Stanton  number.  Experiments  performed  with  the  pressure  fairing  in  place  produce 
results  which  agreed  with  the  previous  correlation  at  zero  pressure  gradient  but  deviate 
from  it  for  the  AT  >  0  cases.  Fig.  3.14  shows  a  comparison  of  results  taken  at  S'  =  0, 
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with  and  without  the  pressure  fairing  in  place.  A  curve  fit  to  aU  of  the  data  at  K  = 
0  had  the  form 


Re  -9.3x10’ — 4-— 

- =  1  -  0.50(1 -e 


(3.7) 


This  is  very  similar  to  the  curve  whidi  was  fitted  to  results  with  the  pressure  fairing 
absent  (Equation  3.3).  Uncertainty  analysis  suggested  that  the  accuracy  to  which 
Stanton  number  and  Re^VRe^^  could  be  measured  was  dependent  on  the  flow  speed 
and  the  power  required  for  a  given  crystal  to  show  colour.  Hie  error  bars  shown  in 
Fig.  3.14  are  for  measurements  of  transition  end  for  each  LC  at  a  flow  speed  of  24.2 
m/s.  They  are  characteristic  of  typical  uncertainties.  Also  shown  is  a  95%  confidence 
interval  for  the  correlation  (Eqn.  3.7).  Details  of  the  uncertainty  analysis  are  given 
in  Appendix  B. 


In  Figure  3.15,  all  results  taken  with  the  pressure  fairing  in  place  (AT  ^  0)  are 
superinqxised  on  the  results  taken  with  the  fairing  absent.  Note  that  all  zero  pressure 
gradient  results,  whether  taken  with  the  fairing  in  place  or  not,  are  given  the  symbol 
in  this  figure.  Under  a  favourable  pressure  gradient,  heating  moves  transition 
forward  a  proportionally  greater  distance  than  an  equivalent  heating  with  zero  pressure 
gradient.  As  AT  is  iiKteased  the  influence  is  larger. 


Figure  3.15  suggests,  however,  that  the  destabilising  influence  of  heating  upon 
the  boundary  layer  is  not  strictly  a  function  of  K,  but  also  varies  with  oncoming  speed 
(lower  oiKoming  qreeds  showing  a  greater  influence  due  to  heating).  Hence,  heating 
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is  shown  to  have  the  greatest  effect  on  fee  normalised  transition  location  at  condition 
Ic  (the  lowest  speed),  although  K  was  greater  for  conditions  2d,  3c  and  3d. 

For  AT  =  0,  aormalised  Reynolds  numbers  at  transiticiin  start  and  end  effectively 
colli^sed  onto  a  single  curve  plotting  against  Stanton  number.  However,  the  results 
do  show  a  trend  towards  a  smaller  proportional  change  in  Re^Re^,^  for  a  given 
Stanton  number  at  transition  end  than  at  transition  start.  This  trend  is  accentuated  for 
K>Q. 


The  length  of  transition,  i)  could  not  be  measured  directly  using  this  heat 
transfer  technique.  As  discussed  in  Section  3. 1.1.1,  the  amount  of  heat  required  to 
identify  transition  end  is  greater  than  that  required  to  identify  transition  start. 
Therefore,  the  transition  start  location  and  associated  heating  rate  for  each  of  fee  LC’s 
was  recorded.  Then  for  each  LC,  the  Reynolds  number  at  transition  end  was 
determined  for  fee  sartK  heating  rate  as  that  used  to  determine  the  location  of 
transition  start  by  interpolation  from  plots  of  transition  end  Reynolds  number  versus 
Stanton  number,  such  as  those  in  Figures  3.9-3.13.  The  transition  length  at  the 
designated  heating  rate  was  then  determined. 


Transition  length  was  found  to  increase  with  increasing  K  as  shown  in  Fig 
3.16.  This  is  consistent  with  the  findings  of  otiiers,  such  as  GAUGLER  (1984). 
Heating  had  no  discetiuble  effect  upon  t|  in  these  experiments. 
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Transition  lengths,  noimalised  with  the  distance  to  die  start  of  transition,  are 
plotted  as  a  function  of  the  heating  rate  in  Fig.  3.17.  Because  of  the  technique  used 
to  find  transition  length,  the  uncertainty  in  q/X,  was  found  to  be  quite  large  (a 
maximum  of  ±  24.0%).  .^>p«idix  B  describes  the  uncertainty  analysis. 

These  results  indicate  that  there  is  little  variation  in  the  normalised  transition 
length  with  either  heating  rate  or  pressure  gradient.  However,  the  results  ate  for  quite 
small  values  of  K  (maxinuim  2.7  x  10^).  Influences  nuty  be  larger  at  higher  pressure 
gradients. 


3.1,2  Carbon-Fibre  Heater 

Figures  3.18  and  3.19  show  heat  transfer  results  from  the  carbon  plate  at  flow 
speeds  of  34.3  and  7.9  m/s  respectively.  The  lines  on  these  plots  represent  the 
following  boundary  layer  heat  transfer  cortelatior.s,  corrected  for  unheated  starting 
length,  from  KAYS  and  CRAWFORD  (1980)  for  flow  over  a  flat  plate  with  a 
constant  heat  flux  and  zero  pressure  gradient: 


0.453 


{laminar)  , 


(3J1) 


0.030ite~°^f>--°* 


{turbulent) 


(33») 


The  unheated  starting  length,  ^  was  measured  as  0.115  m.  Although  INCROPERA 
and  DeWITT  (1985)  note  that  these  expressions  may  be  in  error  by  as  much  as  25% 
due  to  variations  in  surface  roughness  and  free-stream  turbuloKe,  in  tests  similar  to 
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these  SOHN  et  al.  (1989)  and  KIM  et  al.  (1989)  have  had  results  which  agree  with 
these  conelatioiis  quite  well.  The  coirections  for  unheated  starting  lengdi  may  be 
expected  to  break  down  near  the  origin  of  the  thermal  boundary  layer  (see  KAYS  and 
CRAWFORD,  and  TAYLOR  et  al.  (1990)). 

The  level  of  Stanton  numbers  and  their  distributions  suggest  that  the  bouiKlaty 
layer  flow  over  the  heated  section  of  the  plate  was  fully  turbulent  for  the  tests 
presented  in  Figs.  3.18  and  3.19.  Tests  using  the  aluminium  heater,  however,  showed 
the  boundary  layer  to  be  transitional  and  fully  laminar  at  firee-stream  speeds  of  34.3 
and  7.9  m/s  respectively.  Experiments  using  shear  sensitive  liquid  crystal  were  used 
to  investigate  the  state  of  the  bouiulary  layer.  The  strip  of  shear  sensitive  liquid 
crystal  was  60  mm  wide  and  began  4S  mm  downstream  of  the  leading  edge.  It  was 
iq^lied  down  the  centre-line  of  the  carbon  heater  using  the  technique  discussed  in 
Section  2.3.1.  Thb  tests  indicated  diat  the  boundary  layer  may  have  been  disturbed 
by  the  30  pm  step  where  the  perspex  plate  joined  the  catbcm  sattq>le.  The  sequence 
of  photognqdts  presoued  in  Fig.  3.20  shows  the  shear  sensitive  LC  Sowing  Srst  at  the 
juncture  between  the  carbon  and  perspex  plates.  However,  the  results  also  show  the 
shear  sensitive  LC  Sowing  ahead  of  the  junctkm  shortly  afterward.  If  the  boundary 
layer  had  been  laminar  ahead  of  the  junction,  the  shear  stress  should  not  have  been 
high  eiKNigh  to  canae  the  LC  to  Sow. 

Tests  were  also  performed  widi  a  roughness  trip  near  the  leading  edge  in  order 
to  ensure  that  the  boundary  layer  was  folly  turbulem.  The  trq)  consisted  of  a  10  mm 
by  240  mm  str^  of  P40  emery  cloth  which  was  secured  4  mm  downstream  of  the 


leading  edge  using  double-sided  tape.  Figures  3.21  and  3.22  show  heat  transfer  results 
at  flow  speeds  of  28.9  and  7.9  nt/s  respectively.  Measurentents  of  the  coefScimt  of 
heat  transfer  were  generally  insensitive  to  the  presence  of  the  boundaiy  layer  trip. 
This  indicates  that  the  unttipped  boundary  layer  was  also  turbulent  on  the  measure¬ 
ment  surface. 

This  turbtilent  boundaiy  layer  was  attributed  to  blockage  on  the  underside  of 
the  carbon  plate.  At  a  distance  90  mm  downstream  of  the  leading  edge,  the  electrical 
contacts  protiuded  28  mm  into  the  feee-stieam.  As  a  result  of  this  obstruction,  the 
plate  was  effectively  placed  at  incidence  to  the  firee-stieam.  Thus,  the  boundaiy  layer 
separated  at  the  leading  edge,  reattaching  as  a  turbulent  boundaiy  layer. 

Although  the  boundary  layer  was  not  transitional  for  these  tests,  the  results 
suggest  that  the  caibon  heater  may  provide  a  suitable  means  of  obtaining  transition 
locations  and  heat  transfer  measurements.  The  caibon  heater  was  cleady  capable  of 
sustaining  tenqierature  gradienrs  in  the  Itmgitudinal  direction,  as  shown  by  the 
contours  of  constant  tenq)eratuie  which  were  present  during  tests  taken  on  the  fiilly 
turbulent  boundaiy  layer  (see  Fig.  3.23).  Fuithermore,  the  levels  of  Stanton  numbers 
indicated  by  each  of  the  duee  liquid  crystals  are  in  good  accord  and  show  acceptaUe 
agreement  with  standard  cotrelatioas  for  folly  turbulent  flows. 

In  order  to  obtain  transition  on  the  caibon  heater,  the  disnqxion  to  the  flow  on 
the  undeiside  of  foe  plate  was  minimised  by  placing  a  streamlined  jacket  over  the 
electrical  contacts.  The  [date  was  also  elevated  45  mm  at  foe  leading  edge  and  37  mm 
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at  the  trailing  edge,  effectively  creating  an  angle  of  attack  of  -1°.  As  a  tesult, 
transition  was  achieved  on  the  heating  element  (see  Figure  3.24).  The  acceleration 
parameter,  AT,  was  estimated  to  be  2.9  x  lOr*  based  on  measurements  from  diree  static 
pressure  ports  which  were  in  the  perspex  plate  downstream  of  the  carbon  insert. 

Results  obtained  using  this  configuration  are  shown  in  Fig.  3.2S  for  an 
onctmiing  flow  speed  of  28.3  m/s.  Also  presented  is  a  STANS  predictitxi  with 
transition  stan  set  upstream  of  the  heated  section.  The  levels  of  Stanton  number 
indicated  by  each  of  the  three  crystals  are  in  good  agreement.  Calculated  Stanton 
numbers  at  transition  end  are  -  15%  hi^  in  conqrarison  to  the  prediction  but 
iqrproach  the  prediction  in  the  turbulent  region.  Results  from  a  number  of  other 
studies  (e.g.  BLAIR  and  WERLE  (1981),  BLAIR  (1982),  SOHN  et  al.  (1989)  and 
TAYLOR  et  al.  (1989))  show  a  similar  discrepancy  between  predicted  and  experimen¬ 
tal  Stanton  numbers  at  transition  end.  However,  erqretimrats  on  the  aluminium  heater 
displayed  no  such  trend.  Ihis  result  is,  therefore,  most  likely  a  consequence  of  the 
unheated  starting  l«igth.  TAYLOR  et  al.  (1990)  have  found  that  neither  integral  nor 
numerical  boundary  layer  heat  transfer  solutions  accurately  predict  Stanton  numbers 
in  the  heated  region  immediately  downstream  of  an  unbeated  starting  length.  They 
report  that, 

...  a  step  in  wall  leuiperature  has  a  large  effect  on  tire  Stanton  number 
in  the  heated  regitm  near  the  step.  But  as  the  theimal  boundary  layer 
develops,  tite  Stanton  numbers  qjproacb  the  results  for  the  constant 
wall  temperature  boundary  layers. 

Further  confirmation  that  the  boundaiy  layer  was  transitional  was  obtained  by 
placing  a  spherical  excrescence  2  mm  in  diameter  on  dw  surface  of  the  plate  at  a 
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location  120  mm  from  Ok  leading  edge.  Figure  3.26  shows  a  turbulent  wedge  whidi 
formed  behind  Ok  surfoce  roughiKSS  in  the  otherwise  laminar  boundary  layer.  The 
wedge  extended  downstream,  eventually  merging  into  the  turbulent  boundary  layer  at 
the  position  of  natural  transition. 

The  tiitK  requited  for  the  carbon  phue  to  teach  dKrmal  stability  was  estimated 
using  the  method  detailed  in  Section  3.1. 1.1.  This  is  particularly  relevant  to  flight 
Ksting  of  the  NLF  nacelle  as  experimental  titiK  will  be  limited.  Calculated  heat 
transfer  distributions  for  fully  turbulent  boundary  layers  at  speeds  of  28.9  atxl  7.9 
m/sec  are  shown  in  Figs.  3.27  and  3.28  at  titiKS  ranging  from  otK  to  ten  minutes. 
These  results  suggest  that  -  four  minutes  are  required  for  the  carbon  plate  to  reach  full 
thermal  equilibrium,  but  that  quite  accurate  results  can  be  obtairted  after  as  little  as 
two  minutes.  Similar  testing  on  the  aluminium  heater  (Section  3. 1.1.1)  indicated  that 
this  period  may  be  reduced  when  taking  data  in  regions  in  which  Ae  temperature 
gradient  is  high  (i.e.  in  transitional  regions). 

3.U  PreUmmary  Results  fivm  o  Curved  Carbon-Fibre  Heater 

Figure  3.29  shows  a  carbon-fibre  heater  set  into  a  curved  carbon  substrate 
repiesentative  of  a  section  of  the  NLF  nacelle.  The  carbon  structure,  the  heater  and 
electrical  contacts  were  all  manufiKtured  by  Hutel-Dubois  ltd.  to  the  specifications  of 
Rolls-Royce  pic.  and  Oxford  University  in  the  sattK  manner  as  discussed  in  Section 
3.1.2.  The  heater  is  100  mm  x  62S  mm  and  begins  70  nun  from  Ok  leading  edge  of 
the  nacelle  section.  The  model  expands  slightly  in  width  with  distance  from  the 
leading  edge,  measuring  283  mm  (nominally)  x  930  mm.  The  carbon  substrate  curls 


3.15 


under  to  create  a  rounded  leading  edge  in  accordance  with  the  design  of  the  NLF 
nacelle.  Except  for  the  heater  (which  remains  black),  the  nacelle  sectitm  is  white  (as 
required  for  the  thermal  imaging  transition  detection  technique  to  be  performed  by 
DLR).  The  white  coating  is  feathered  out  at  the  edges  of  the  heater  such  that  no 
detectable  step  exists  between  the  heater  and  the  surrounding  carbon  structure. 
Resistance  temperature  detectors  (RTD’s)  25mm  x  70  mm  were  mounted  along  the 
underside  of  the  heater  at  100  mm  increments.  For  these  tests,  the  underside  of  die 
heater  was  shielded  from  the  free-stream  flow  by  a  loose  plastic  covering,  thus 
providing  an  insulating  regirni  of  dead  air.  The  model  was  mounted  vettically  in  the 
tunnel  such  that  it  was  flush  with  the  upper  and  lower  surfaces  of  the  tunnel  test 
section  (refer  to  Fig.  3.29).  The  angle  of  attack  of  the  nacelle  section  was  adjusted 
by  rotating  it  about  a  mounting  screw  attached  midway  along  its  length  and  extended 
through  the  tunnel  floor. 

The  heater  conformed  well  to  the  double  curvature  of  the  nacelle  section.  An 
imperfect  contact  at  the  upstream  end  of  the  heater  led  to  a  hot-spot  there.  This  did 
not  preclude  a  qualitative  test  of  the  heater  performance.  Figure  3.30  shows  the 
nacelle  section  with  a  fully  laminar  boundary  layer  along  the  length  of  the  heater.  A 
turbulent  wedge  formed  downstream  of  an  excrescence  attached  to  the  surfoce  of  the 
heater,  thereby  confirming  the  lamiruu  luture  of  the  boundary  layer.  The  RTD’s 
attached  to  die  backside  of  the  heater  acted  as  small  heat  sinks  (as  shown  in  Figure 
3.30).  However,  their  presence  is  not  thought  to  offer  a  significant  inqiedimern  to  heat 
transfer  tests. 
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By  vaiying  the  angle  of  attack  of  the  nacelle  section,  it  was  also  possible  to 
cause  the  boundary  layer  to  undergo  transition  along  the  test  surface.  Hgure  3.31 
shows  the  location  of  transition  start  as  indicated  by  an  LC.  The  hot-spot  due  to  die 
inqierfect  contact  is  s^^iarent  at  the  upstream  electrical  contact.  No  heating  was 
detected  at  the  downstream  contact. 

These  qualitative  results  represent  only  preliminary  investigations.  Quantitative 
measurements  of  transition  and  heat  transfer  on  the  nacelle  section  are  cuirently  under 
way  at  Oxford  University. 

3.2  Shear  Sensitive  Liquid  Crystal 

3.2.1  Natural  Transition 

Shear  sensitive  liquid  ciystal  was  used  qualitatively  to  deteimine  the  end  of 
transition  by  noting  the  location  at  which  the  liquid  crystal  film  first  changed  colour. 
A  direct  comparison  of  transition  end  location  as  determined  Ctom  the  heat  transfer 
experiments  and  transition  location  as  determined  firom  the  shear  stress  experiments 
was  obtained  by  placing  a  60  mm  str^  of  shear  sensitive  LC  alongside  die  aluminium 
heating  str^.  This  necessarily  placed  the  shear  sensitive  LC  towards  the  side  of  the 
plate  where  experiments  indicated  that  cross-flow  was  present  in  die  boundary  layer. 
The  location  of  transition  end  as  determined  by  extrapolating  the  heat  transfer  results 
to  adiabatic  conditions,  agreed  with  the  location  of  transition  eiMl  determined  from  the 
shear  experiments  to  within  15%  over  the  range  of  speeds  at  which  the  shear  LC  could 
be  used. 
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More  detailed  tests  using  the  shear  sensitive  LC  were  peifoimed  on  a  plate 
with  a  slightly  different  leading  edge.  The  tests  were  conducted  at  die  high  end  of  the 
tunnel's  speed  range,  where  the  surface  shear  stress  was  relatively  high  (but  less  than 
5  Pa).  At  lower  tunnel  speeds  the  shear  stress  indication  was  influenced  to  a  greater 
extent  by  non-unifoimities  in  the  layer  of  shear  sensitive  LC.  For  levels  of  shear 
stress  <  2  Pa  the  LC  film  changed  colour  in  patches,  thereby  rendering  the  method 
unreliable  for  quantitative  measurements.  Reliable  results  were  obtained  at  36.5  and 
33.4  m/s.  BONNETT  et  al.  (1989)  have  repotted  that  it  may  be  possible  to  obtain 
results  at  lower  shear  stresses  by  reducing  the  viscosity  of  the  LC. 

A  sequence  of  photogr^hs  showing  the  LC  layer  at  different  times  at  the 
highest  flow  speed  is  presented  in  Figure  3.32.  The  first  photograph  after  the  flow 
starts,  r  =  19  sec.,  shows  the  location  of  the  end  of  transition  (highest  shear  stress)  at 
about  230  nun  finm  the  leading  edge  of  the  plate.  At  greater  times,  the  location  at 
which  the  layer  starts  to  flow  moves  downstream  rrqiidly  and  upstream  slowly.  This 
is  associated  with  the  small  gradient  in  shear  stress  in  the  turbulent  boundary  layer 
conqiated  with  that  in  the  transitional  region. 

BONNHTT  et  al.  show  that  the  type  of  cholesteric  liquid  crystal  used  has  a 
calibration  of  the  form 

ln(T)  =  A +BIn(f)  .  (3-10) 

where  t  is  the  surface  shear  stress,  t  is  the  tuiK  for  the  LC  to  show  colour  and  A  and 
B  are  constants  found  by  calibradrm.  An  attenqrt  to  calibrate  this  LC  has  been  made 
by  fitting  the  experimental  data  in  the  turbulent  region  for  the  highest  speed  runs  to 
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a  prediction  from  STANS.  The  result  of  such  a  procedure  is  shown  in  Fig.  3.33, 
where  the  values  determined  for  A  and  B  are  1.9  and  -0.18  respectively.  Also  shown 
are  the  experimental  results  from  the  lower  speed  runs  plotted  using  that  calibration 
and  compared  with  anodier  prediction.  In  both  cases,  the  experimentally  determined 
transition  location  has  been  used  for  the  prediction.  The  results  are  encouraging,  but 
further  tests  over  a  wider  range  of  shear  stress  are  required  for  validation. 

At  present  the  teduiique  requires  the  observer  to  assess  the  locations  at  which 
the  layer  of  LC  has  started  to  move  at  a  given  point  in  dme.  This  requires  some 
subjective  judgement  if  the  layer  is  not  uniform.  However,  the  form  of  the  calibration 
(Eqn.  3.10),  is  such  that  a  30%  error  in  t  typically  leads  to  only  a  6%  error  in  r. 
Experiments  also  indicate  that  the  time  required  for  the  LC  to  show  colour  is 
influenced  by  temperature.  This  is  associated  with  the  tenqwiatute  dependence  of  the 
viscosity  of  the  mixture. 

The  influence  of  temperature  is  qrparem  in  Figs.  3.34  and  3.3S,  which  show 
the  repeatability  of  measurements  using  shear  sensitive  LC.  The  same  calibration 
constants  were  used  for  all  tests  presented  regardless  of  the  operating  tenqrerature. 
Repeated  tests  at  a  given  speed  consistently  indicated  the  same  location  for  transition 
end. 


Some  tests  were  performed  by  WALTON  (1991)  with  a  fully  developed  turbu¬ 
lent  flow  in  a  duct  of  cross-section  3  mm  x  7S  mm  in  an  attempt  to  independently 
calibrate  the  same  cholesteric  liquid  crystal  of  the  same  thickness.  Preliminary  results 
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from  these  tests  indicate  that  the  form  of  the  calibration  is  as  in  Eqn  3.10.  Slightly 
different  values  of  A  and  B  were  indicated  (1.75  and  -0.175),  but  die  values  of  shear 
stress  obtained  with  this  calibration  were  within  15%  of  those  in  Fig.  3.33. 

The  method  of  using  shear  sensitive  liquid  crystal  will  not  be  suitable  for  use 
in  regions  with  large  negative  gradients  of  shear  stress  in  the  direction  of  flow.  The 
LC  flows  much  more  rapidly  once  it  has  entered  the  Grandjean  texture,  which  it  wiU 
do  at  earlier  times  in  regions  of  higher  shear  stress.  The  LC  wiU,  therefore,  flow 
downstream  over  parts  of  the  LC  layer  where  shear  stress  is  low  and  the  layer  remains 
in  the  focal  conic  texture.  Thus,  results  will  be  masked  in  areas  of  low  x. 

3.2,2  Induced  Transition 

In  order  to  insure  that  the  shear  sensitive  liquid  crystal  was  performing  as 
expected,  small  excrescences  were  attached  to  the  plate  55  mm  from  the  leading  edge 
to  produce  turbulent  wedges  in  an  otherwise  laminar  boundary  layer.  Sequences  of 
photographs  showing  the  development  of  the  wedges  at  two  flow  speeds  are  presented 
in  Figures  3.36  and  3.37.  For  these  tests,  the  flow  was  started  without  any 
excrescences  and  the  LC  was  allowed  to  change  colour  in  the  transition  region.  The 
flow  was  then  stopped.  The  excrescence  was  attached  to  the  surface  of  the  plate  and 
the  tunnel  was  restarted.  For  the  higlMr  speed  test,  the  excrescence  height  was  0.3 
mm.  This  was  found  to  be  insuffleient  to  generate  a  wedge  at  the  lower  speed,  and 
the  height  of  the  excrescence  was  increased  to  0.5  mm. 

These  tests,  particularly  that  at  the  lower  speed,  show  the  variation  of  shear 
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stress  within  the  wedge.  Contours  of  shear  stress  are  presented  in  Figure  3.38.  As 
would  be  expected,  the  shear  is  higher  close  to  the  excrescence  where  the  hydrostatic 
boundary  layer  is  thinner.  The  semi-angle  of  the  wedge  increases  with  distance 
downstream  as  has  been  observed  previously  for  small  excrescences  by  SCHUB  AUER 
AND  KLEBANOFF  (19SS).  The  magnitude  of  the  semi-angle  (approximately  6°) 
suggests  that  the  shear  LC  has  changed  colour  only  in  the  fully  turbulent  region  of  the 
wedge  where  the  shear  stress  will  be  hipest  (see  Fig.  3.39).  The  results  also  indicate 
that  the  shear  stress  is  higher  towards  the  edges  of  the  wedge  than  near  the  centre. 
Considering  the  variation  of  shear  stress  along  a  span-wise  line  it  is  noted  that  the 
level  of  shear  stress  rises  from  the  laminar  value  outside  of  the  wedge  and  increases 
through  an  intermittent  region  at  tire  edge  of  the  wedge  (although  the  fuU  extent  of 
this  rise  has  not  been  ctptuted  in  the  present  experinrent).  The  shear  stress  attains  a 
maximum  value  at  the  edge  of  the  folly  turbulent  region  of  the  wedge  but  decreases 
toward  the  centre. 

These  results  are  in  general  agteenrent  with  those  of  CLARK  and  JONES 
(1991)  who  obtairred  contours  of  heat  transfer  in  a  turbulent  wedge  using  the  steady- 
state  heat  transfer  technique  and  temperature  sensitive  liquid  crystals.  They  have 
noted  that  the  level  of  heat  transfer  rises  from  the  laminar  value  outside  the  wedge  and 
increases  through  an  intermittent  region  at  the  edge  of  the  wedge  (see  Fig.  3.40). 
However,  CLARK  and  JONES  have  not  taken  data  for  semi-angles  of  <  3°  (the  fully 
turbulent  portion  of  the  wedge).  Therefore,  they  have  not  reported  the  decreasing  rate 
of  heat  transfer  along  a  span-wise  litre  in  dre  frilly  turbulent  section  of  the  wedge, 
which  was  indicated  by  die  shear  crystal. 
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Longitudinal  streaks  are  also  apparent  in  the  LC  film,  especially  early  in  the 
sequence  of  photographs  at  the  lower  speed.  These  ^>pear  to  be  associated  with  the 
presence  of  stream-wise  vortices.  THOMSON  (1988)  has  reported  that  stream-wise 
vortices  wrtq)  around  the  protuberance  and  trail  downstream.  The  location  in  planfoim 
of  six  vortex  legs  can  be  distinguished  in  the  LC  film  (first  photograph  of  Fig  3.37). 
Similar  flow  visualisation  has  been  achieved  using  sublimation  of  aceniq>hthene  on  a 
flat  plate,  as  shown  in  Fig.  3.41.  THOMSON  suggests  that  the  induced  velocities 
associated  with  these  vortices  is  responsible  for  flow  breakdown  to  turbulence.  Hence, 
turbulence  appears  to  occur  first  directly  downstream  of  the  surface  roughness  in  the 
vicinity  of  the  innermost  vortex  legs  which  are  the  strongest.  Transition  occtus  in  the 
proximity  of  the  other  vortex  legs  further  downstream.  The  LC  film’s  capability  to 
detect  stream-wise  vortices  suggests  that  the  shear  sensitive  LC  is  quite  sensitive  to 
variations  in  surface  shear  stress. 

These  results  show  the  possible  usefulness  of  shear  sensitive  liquid  crystals  for 
fiill  coverage  measurement  of  surface  shear  stress  in  transitional  boundary  layer  flows. 
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Chapter  4:  Conclusions  and  Recommendations 


The  present  study  is  primarily  diteaed  at  developing  a  steady  state  heat 
transfer  technique  employing  temperature  sensitive  liquid  crystals  for  transition 
detection.  Conclusions  related  to  the  technique  itself  are  presented  here,  follovkred  by 
those  related  to  the  development  of  a  constant  heat  flux  heater  constmcted  of 
unidirectional  carbon-fibres.  Finally,  conclusions  are  presented  from  experiments 
employing  liquid  crystal  to  measure  surface  shear  stress. 

4.1  Heat  Tranter  Measurements 

Tests  performed  on  a  flat  plate  in  incompressible  flow  suggest  that  the  steady 
state  heat  transfer  technique  employing  enc^ulated  temperature  sensitive  liquid 
crystals  provides  a  suitable  means  for  obtaining  high  resolution  heat  transfer 
measurements  and  detecting  boimdary  layer  transition.  Self-consistent  measurements 
of  heat  transfer  are  obtained  from  a  multiple  LC  layer  and  agree  well  with  standard 
correlations  in  the  laminar  boundary  layer.  Results  fall  below  the  correlations  in  the 
turbulent  portion  of  the  boundary  layer  (particulariy  near  transition  end)  by  a 
maximum  of  13%,  but  qtproach  the  correlations  further  downstream. 

Measured  coefBcients  of  heat  transfer  in  the  laminar  and  turbulent  regions  of 
the  boundary  layer  were  found  to  be  unaffected  by  the  level  of  heat  transfer  from  the 
surface.  However,  results  clearly  indicate  that  heating  the  surface  significantly  hastens 
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transition  and  must,  therefore,  be  considered  when  inferring  the  point  of  transition  for 
an  unheated  surface.  By  using  three  liquid  crystals,  each  with  a  different  colour-play 
temperature,  it  is  possible  to  plot  the  movement  of  transition  as  a  function  of  heating 
for  a  given  condition.  An  exponential  extrapolation  to  zero  heating  then  provides  a 
measure  of  the  adiabatic  transition  Reynolds  number.  Plotted  against  a  normalised 
transition  Reynolds  number,  a  Stanton  number  of  the  form  q/(puC^o)  correlates  the 
movement  of  transition  with  heating.  Uncertainties  in  measurements  of  Stanton 
numbers  and  normalised  transition  Reynolds  numbers  are  a  maximum  of  9.5%  and 
6.9%  respectively  in  the  present  tests. 

The  destabilising  influence  of  heating  on  the  flat  plate  boundary  layer  is 
counteracted  by  the  stabilising  influence  of  a  favourable  pressure  gradient.  Under 
adiabatic  conditions,  the  transition  Reynolds  number  increases  rapidly  with  favourable 
pressure  gradient  and  a  single  ’acceleration’  parameter,  K,  is  found  to  correlate 
Reynolds  number  at  transition  start  and  end.  With  favourable  pressure  gradients, 
heating  has  a  greater  influence  on  transition  location  dian  in  zero  pressure  gradient. 
As  K  increases,  the  infltience  of  heating  increases.  The  present  results  suggest  that 
a  single  acceleration  parameter  alone  does  not  correlate  the  movement  of  transition 
with  heating.  With  zero  pressure  gradient,  normalised  Reyrtolds  numbers  at  transition 
end  effectively  collapse  onto  a  single  curve  plotting  against  Stanton  number.  With 
favourable  pressure  gradient,  results  show  a  trend  towards  a  smaller  change  in 
normalised  Reynolds  numbers  with  heating  at  transition  end  than  at  transition  start. 

Using  the  present  technique  with  luirrow-band  liquid  crystals,  transition  length 
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can  not  be  measured  diiecdy,  but  can  be  extracted  indirectly.  Transition  length  is 
found  to  increase  with  increasing  K,  and  heating  has  no  discernible  effect  on  transition 
length.  Normalised  transition  length  varies  little  with  heating  or  the  small  pressure 
gradients  of  the  present  experiments.  The  method  of  obtaining  transition  length 
introduces  significant  unceitainties  (a  maximum  of  ±  24%  in  the  noimalised  transition 
length).  Hieses  unceitainties  could  be  reduced  significantly  by  using  a  wide-band  LC 
for  transition  detection  (as  discussed  below). 

4.1.1  Future  Work 

Narrow-band  temperature  sensitive  liquid  crystals  have  been  used  in  the  present 
heat  transfer  experiments.  These  enable  surface  tenqseratuies  to  be  accurately 
measured  which,  in  turn,  allows  heat  transfer  coefficients  to  be  found.  However,  this 
work  is  primarily  a  background  study  into  using  the  steady  state  heat  transfer 
technique  to  measure  transition.  This  requires  only  that  the  positions  of  minimum  and 
maximum  surface  temperatures  be  located  for  a  number  of  heating  rates  in  order  to 
infer  the  adiabatic  transition  location.  This  may  be  more  efficiently  accomplished  by 
using  a  wide-band  LC  which  b  capable  of  indicating  the  location  of  transition  start 
and  end  at  a  single  heating  rate  and  (unlike  narrow-band  LC’s)  will  not  require  the 
heating  rate  to  be  adjusted  so  as  to  position  a  particular  ist^hermal  band  at  the  location 
of  minimum  or  maximum  surface  tenqierature. 

A  proposed  technique  b  to  use  a  wide-band  LC  (to  locate  transition) 
accompanied  by  a  narrow-band  LC  (to  obtain  accurate  values  of  heat  transfer 
coefficients).  Accurate  surface  temperature  measurements  can  be  supplemented  by 
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resistance  temperature  detectors  (RTD’s)  but,  of  course,  these  will  only  provide 
measurements  at  discrete  locations.  Tests  with  combined  wide-  and  narrow-band  LC 
coatings  ate  under  way  at  O.U.E.L. 

4.2  Vnidinctional  Carbon-Fibre  Heater 

A  constant  heat  flux  heater  constructed  from  unidirectional  caibon  fibres  was 
also  tested  in  a  wind  tunnel  is  proposed  for  use  in  in-flight  transition  detection 
experiments.  The  heater  is  to  be  integrated  into  the  load-bearing  structure  of  an 
aircraft  engine  nacelle.  To  the  audior’s  knowledge  this  represents  the  first  attempt, 
to  locate  transition  and  make  heat  transfer  measurements  employing  part  of  the  aircraft 
stracture  as  a  heating  element.  Such  an  experimental  technique  may  prove  very  useful 
in  a  variety  of  future  flight  test  plications. 

The  unidirectional  caibon  fibres  produce  an  extremely  uniform  surface  heat 
flux.  There  is  no  detectable  step  between  the  heater  and  the  surrounding  carbon 
structure.  The  heater  is  ctqiable  of  sustaining  the  lateral  and  longitudinal  temperature 
gradients  necessary  for  detecting  transition  and  making  heat  transfer  measurements 
using  the  present  technique.  Measurements  of  heat  transfer  in  a  tuibulent  boundary 
layer  obtained  using  the  carbon  beater  agree  weU  with  standard  coneladons.  Limited 
tests  in  transitional  boundaiy  layers  found  measuiements  of  heat  transfer  near  the  end 
of  transition  to  be  a  maximum  of  -  15%  higher  than  STANS  predictions.  This  is 
attributed  to  boundaiy  layer  conditions.  Further  downstream  (in  the  well-developed 
portion  of  the  turbulent  boundary  layer),  heat  transfer  measuiements  ap|Hoach 
predictions.  The  caibon  plate  requires  approximately  four  minutes  to  readi  full 
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thermal  equilibrium  in  the  present  tests,  although  quite  accurate  results  can  be  obtained 
after  as  little  as  two  minutes. 

Good  electrical  contact  to  the  carlxMi  heater  is  essential,  as  poor  contact  results 
in  non-uniform  heating  of  the  test  surface  and  contact  heating.  Suitable  electrical 
contacts  have  been  devised  and  successfully  tested.  Resistance  in  the  exposed  carbon 
fibres  leading  to  the  electrical  contacts  may  be  effectively  eliminated  by  soaking  the 
fibres  in  a  mixture  of  silver-loaded  conducting  paint  and  xylene,  then  allowing  the 
xylene  to  evtqwrate. 

The  thermal  conductance  of  the  carbon-fibre  matrix  is  estimated  to  be  6.7 
Wm  '  K  ‘.  Due  to  this  relatively  high  thermal  conductivity,  adequate  thermal 
insulation  on  the  underside  of  the  carbon  beater  (to  include  the  contacts  and  the  area 
surrounding  them)  is  very  inqwrtant.  Conservative  estimates  suggest  that  lateral  and 
longitudinal  conduction  may  amount  to  a  maximum  of  2-8%  and  13%  respectively  for 
these  tests. 

Preliminary  tests  on  a  carbon-fibre  heater  set  into  a  curved  section 
representative  of  that  proposed  for  the  NLF  nacelle  indicate  that  transition  can  be 
detected  using  teiiqwratute-sensitive  liquid  crystals. 

43  Sutface  Shear  Stress  Measurememts 


In  related  experiments,  shear  sensitive  liquid  crystal  is  used  for  the  first  time 


to  obtain  quantitative  measuiements  of  surface  shear  stress  in  a  transitional  boundary 
layer.  The  technu]ue  relies  upon  calibrating  the  time  required  for  an  LC  layer  to 
change  from  the  focal-conic  texture  to  the  Grandjean  texture  when  subjected  to  a 
shearing  force. 

There  is  a  lower  limit  of  shear  stress  (2  Pa)  below  which  reliable  results  can 
not  be  obtained  using  the  present  materials.  Therefore,  measurements  could  not  be 
obtained  in  the  laminar  portion  of  the  boundary  layer  with  the  flow  speeds  available 
in  the  present  tunnel. 

The  shear  sensitive  IX!  calibration,  obtained  in  the  test  rig  using  a  STANS 
prediction,  was  found  to  lead  to  results  at  another  test  condition  in  good  agreement 
with  another  STANS  prediction.  Independent  calibrations  produced  results  within  1S% 
of  those  which  used  the  calibration  firom  the  test  tig. 

Repeated  tuns  at  similar  conditions  provide  consistent  indications  of  transition 
end  location.  Transition  end  Reynolds  numbeis  as  given  by  shear  sensitive  LC  agree 
with  the  (adiabatic)  transition  end  Reynolds  numbers  obtained  usirtg  the  heat  transfer 
technique  to  within  1S%  when  both  tests  ate  performed  on  the  same  plate. 

Shear  sensitive  LC  was  also  used  to  obtain  contours  of  shear  stress  within  a 
turbulent  wedge  induced  by  a  surface  excrescence.  Results  were  obtained  in  the  fully 
turbulent  region  of  the  wedge  atxl  partially  in  the  intemiittent  regions  on  either  side 
of  the  wedge.  In  the  span-wise  direction  the  level  of  shear  stress  rises  from  the 
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laminar  value  outside  the  wedge  and  is  shown  to  increase  through  the  intermittent 
region  at  the  edge  of  the  wedge  until  reaching  a  maximum  value  at  die  start  of  the 
fully  turbulent  region.  Shear  stress  then  decreases  toward  the  centre  of  the  wedge. 
Stream-wise  voitices  downstream  of  the  excrescence  were  also  visible  in  the  LC  film. 

These  results  show  the  potential  usefulness  of  shear  sensitive  LC’s  for  full 
coverage  measurement  of  surface  shear  stress  in  transitional  boundary  layers.  While 
results  are  encouraging,  fiitdier  tests  over  a  wider  range  of  shear  stress  levels  are 
required  for  validation.  Shear  sensitive  LC  will  not  be  suitable  for  use  in  regions  of 
large  negative  gradients  of  shear  stress  in  the  flow  direction. 
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Appendix  A:  Dimensional  Analysis 

Consider  the  flow  of  an  incompressible  fluid  past  a  uniformly  heated  flat  plate 
on  which  a  boundary  layer  forms  and  undergoes  transition  (see  Fig.  A.l).  The 
physical  quantities  which  determine  the  location  of  transition  are:  the  free  stream 
velocity,  u,  the  viscosity,  ft,  the  density  of  the  fluid,  p,  the  gas  ten^rature,  T^,  the 
specific  heat  of  the  fluid  at  constant  pressure,  C^,  a  length  scale,  /,  the  distance  to  the 
start  of  transition,  x,,  the  thermal  conductivity  of  the  fluid,  k.  the  acceleration  of  the 
flow,  du/dx,  and  the  heat  flux  from  the  plate,  q. 

The  fundamental  dimensions  required  to  describe  these  10  physical  variables 
are 

M  =  dimensions  of  mass  , 

L  =  dimensitHis  of  length  , 

T  =  dimensions  of  time  , 

0  s  ditrtensicHis  of  tenqx;rature  , 

Q  =  dimensions  of  heat 

The  physical  variables  and  their  dimensions  ate 

(«]  =  LT' 

M  =  MT'L'  , 

ip]  =  wz.’ 
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IT,]  =  e 

[C^  =  GAf'e  ‘  . 
\n  =  L 
M  =  L 

m  =  QT'L  'e  '  , 
[duJdx]  =  T' 

[4]  =  QT'L-^  . 

The  resulting  indicial  matrix  follows. 


quantity 

M 

L 

T 

0 

Q 

u 

0 

1 

-1 

0 

0 

M 

1 

-1 

-1 

0 

0 

P 

1 

-3 

0 

0 

0 

0 

0 

0 

I 

0 

c. 

-1 

0 

0 

-1 

1 

1 

0 

1 

0 

0 

0 

X, 

0 

1 

0 

0 

0 

k 

0 

-1 

-1 

-1 

1 

du/dx 

0 

0 

-1 

0 

0 

4 

0 

-2 

-1 

0 

1 

Note  that  rows  associated  widi  the  jrfiysical  variables  k,  ft,  q,  I,  and  du/dx  are  linearly 
dependent.  Therefore,  the  rank  of  the  indicial  matrix  is  flve. 

Following  the  Buckingham  pi  theorem  as  presented  in  ANDERSON  (1984), 
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let  N  equal  the  number  of  physical  variables  and  R  equal  the  rank  of  the  indicial 
matrix.  Let  P„  P2,  .  .  .  ,  Pn  represent  the  physical  variables  in  the  relation 


/.(P,.  Pj,  .  .  .  .  P^)=0  .  (A.1) 

Then,  the  physical  relation  in  equation  A.l  may  be  re-expressed  as  a  relation  of  (N  - 
R)  dimensionless  11  products, 

^  (  Hj  ,  Hj  ,  .  .  .  ,  Iljy.jj)  =0  , 

where  each  11  product  is  a  dimensionless  product  of  an  arbitrarily  selected  set  of  R 

physical  variables  plus  one  other  physical  variable.  Let  P„  Pj . P,  be  the 

arbitrarily  selected  set  of  R  physical  variables.  Then 

~  . » 

^2  ^  . » 


= /s (^1  .^2' •  ■  ■ • 

Since  N  =  10  and  R  =  5,  equation  A.l  may  be  expressed  as  a  relation  of  five 
n  products.  Using  p,  «,  x^,  k,  and  Tg  as  repeating  variables  the  following 
dimensionless  quantities  emerge: 


Dt  = 


pax, 
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n,  = 


n4  = 


n* 


ft  ^ 

u  dt 


These  dimensionless  quantities  may  now  be  combined  in  various  manners  in  order  to 
produce  other  more  useful  dimensionless  products.  Reynolds  number  is  already 
represented  by  11,.  Multiplying  n,  by  IIj  produces  a  Prandtl  number.  The  product 
of  l/TIj  and  H,  is  a  unit  Reynolds  number.  A  dimensionless  heating  parameter,  or 
Stanton  number,  is  created  by  multiplying  II4  by  1/11,.  Finally,  a  convenient 
acceleration  parameter,  K,  is  formed  by  combining  FIj  and  FI,.  The  new  n  quantities, 
then,  are 


n.  = 


k 


n. 


pul 


n4  - 


puC. 


T 

p  t 
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Appendix  B:  Vncertamty  Analysis 


An  analysis  of  the  unceitainties  in  measurements  of  transition  length  and 
Reynolds  and  Stanton  numbers  at  transition  start  and  end  was  petfonned.  The 
measurements  which  were  analyzed  were  obtained  using  the  aluminium  beater  with 
zero  pressure  gradient  conditions.  Since  the  resistance  of  the  carbon-fibre  heater  is 
more  unifoim  than  that  of  the  aluminium  beater,  uncertainties  in  results  obtained  using 
the  fonner  are  expected  to  be  smaUer.  The  e:q)erimental  arrangement  was  the  same 
regardless  of  which  heater  was  employed. 

FoUowing  KUNE  and  McOJNTOCK  (1953)  and  BAINES  et  al.  (1991),  let 
an  experimental  result  F  be  a  linear  function  of  n  independent  variables,  4>t, 
each  of  which  is  normally  distributed  and  possesses  the  same  level  of  odds 

F  =  F(4>i,  4)2.  ...«!>,)  (B.1) 

Each  nuiy  be  written  as 

♦i  =  ,  (B  J) 

where  is  the  true  value  and  S^,  is  the  uiKertainty.  The  component  of  uncertainty 
of  F  due  to  an  imceitainty  in  4'.  is  ^tcn 

(8F).  = 

The  expression  for  the  uncertainty  of  the  result  which  gives  the  same  odds  as  each  of 
the  variables  is 
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fiF  = 


(B.4) 


^  i6F)l  *  ...  *  (aF)^] 

Unceitainties  in  measurements  of  the  following  variables  were  ctmsidered: 

(1)  electrical  cunent,  i, 

(2)  total  temperature,  Tg, 

(3)  calibrated  tm^rature  at  which  the  liquid  crystal  shows  colour,  T^, 

(4)  location  of  isothermal  line,  x 

(5)  Ctee-stream  velocity,  u, 

(6)  local  resistance  of  die  heater,  F. 

Uncertainty  in  each  individual  variable  may  derive  from: 

(1)  limitations  in  die  precision  to  whidi  an  instrument  can  be  read, 

(2)  unsteadiness  or  non-uniformities, 

(3)  instrument  calibration  inaccuracies. 

Values  of  the  current  being  fed  to  die  heater  were  read  from  a  Salford  Selectest 
model  ’Super  SO’  analog  ammeter.  Unceitainties  in  measurements  of  i  were  due  to 
mininnmi  resolution  of  die  analog  meter  (±  0.02  anqis),  fluctuations  in  cunem 
provided  by  die  variable  voltage  source  (±  0.03  amps)  and  die  calibrated  accuracy  of 
the  ammeter  (0.3%  of  full  scale  ±  2.5%  of  die  reading).*  A  10  amp  scale  was  used 
for  diese  experiments. 

*Publisbed  accuracy  for  die  Salford  Selectest  is  ±  2.5%  of  full  scale.  However, 
experiments  showed  diat,  when  die  Selectest  was  mounted  horizontally  and  was 
monitoring  current  in  the  same  arrangemem  as  for  the  experimem,  its  measurements 
of  current  agreed  with  measurements  from  a  Tfauilby  1S03-HA  digital  multimeter 
widiin  the  accuracy  of  die  latter  instrument  tfarou^iout  the  range  of  curiait  settings 
enqiloyed  in  these  esperimeias.  Iherefoie,  uncertainty  analysis  was  performed  using 
the  published  accuracy  of  the  1503-HA. 
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Measunments  of  mean  free-stieam  total  temperatuie  were  made  using  an 
aspirating  total  temperature  probe  ctmtaining  a  type  T  dieimocouple.  A  Comaik 
model  6900  digital  diennometer  was  used  to  mraiitor  die  ten^ierature  of  the 
thermocouple.  Estimated  ST^  was  ±  0.2°C.  This  contributed  <  0.3%  uncertainty  to 
measured  Stanton  numbers  and  was,  dierefore,  neglected. 

Uncertainties  regarding  the  tenqierature  at  which  a  specific  colour  band  of  a 
given  LC  shows  colour  were  due  to  the  limited  accuracy  of  the  Cmnaik  model  6900 
themKMneter  (±  0.2‘’C),  uncertainty  in  ascertaining  the  location  of  the  isodietmal  line 
in  relation  to  the  calibrating  thermocouple  (±  0.2°C)  and  drifting  in  die  re^xmse 
temperature  of  the  LC  (±  0.2®C).  Total  ST,  was  ±  0.3®C.  This  typically  contributed 
<  0.1%  uncertainty  to  measurements  of  Stantrai  numbers.  Consequendy,  ST^  was 
neglected. 

The  degree  of  certainty  to  which  the  location  of  an  isothermal  ctmtour  could 
be  measured  varied  depending  <hi  whedier  the  contour  designated  the  location  of 
transition  start  or  end.  Temperature  gradients  are  shaiper  in  the  region  of  minimum 
heat  transfer  than  in  the  region  of  maximum  heat  transfer.  Consequendy,  the  colour 
conttnirs  are  shaiper  near  transidon  start.  Near  x,  die  location  of  the  comour  could  be 
determined  to  ±  3  mm,  while  resolution  was  ±  7  mm  near  Xf. 

The  resistance  of  the  metallised  polymer  was  6.29  £2,  as  measured  by  a 
Thuilby  model  1S03-HA  digital  multimeter.  The  iS03-HA  is  calibrated  to  an 
accuracy  of  0.0S%  f.s.  ±  0.15%  of  the  reading.  Uncertainty  in  the  value  of  the  local 
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resistance  due  to  non-uniformities  was  assumed  to  be  ±  5%  (ref.  Sect  2.2.1). 
Theiefoie,  SR  was  ±  0.36  Q. 

Free-stieam  velocity  was  calculated  from  measurements  of  total  and  static 
pressure  in  the  fiee-stream.  The  measurements  were  made  using  a  Pitot-static  tube 
and  a  manometer  bank  inclined  at  30  deg.  Total  pressure  could  be  read  to  within  ± 
0.5  mm  of  manometer  fluid.  Fluctuations  in  tunnel  qieed  added  uncertainty  to  static 
pressure  measurements,  which  could  be  read  to  within  ±  1  mm  of  manonteter  fluid. 
Total  uncertainty  amounted  to  ±  1.1  mm  of  manometer  fluid. 

In  transitional  boundary  layers,  the  largest  uncertainties  in  measurements  of  Re 
and  St  occur  at  low  speeds  where  die  relative  uncertainty  in  i  and  u  is  maximal.  The 
lowest  speed  at  which  measurements  of  transition  start  and  end  were  avaflable  was 
22.3  m/s  (designated  flow  condition  la  in  Section  3.1. 1.2).  TaUe  B.l  shows 

Table  B  J  Transititxi  locations  and  Reymrids  numbers,  and  their  associated  heating 
rates,  at  223  m/s.  Flat  {date,  zero  pressure  gradient. 


Transition  Start 

Transition  End  | 

X, 

(m) 

Re. 

(xlO^) 

St. 

Xt 

(m) 

Ref 

Stf 

(xlO*) 

23.7X 

crystal 

.210 

3.04 

3.39 

.330 

4.78 

5.91 

29.7*C 

crystal 

.200 

2.87 

5.94 

.290 

4.16 

10.7 

37.1°C 

crystal 

.165 

2.35 

9.23 

.225 

3.20 

17.4 
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'  nteasurements  of  transition  start  and  end  locations  and  Reynolds  numbers,  aitd  tiieir 

’  associated  heating  rates,  for  each  LC  at  tiiis  speed. 

I 

I  B.I.  Vncertaiiay  in  Reynolds  and  Stanton  Numbers  at  Transition  End 

I 

Using  E^.  B.3  and  B.4,  uncertainties  in  Reynolds  and  Stanton  numbers  were 
I  calculated  from  die  uncertainties  associated  with  measurements  of  i,  R.  x  and  u. 

j 

I  Results  are  shown  in  Table  B.2  for  transititm  end,  where  uncertainty  is  expected  to  be 

I 

highest  (since  Sxf  >  Sx,). 

/ 


Table  B,2  Uncertainty  in  Reynolds  and  Stanttm  numbers  at  transition  end  on  a  flat 
plate  with  zero  pressuie  gradient.  Speed  =  24.2  m/s. 


B  JS.  Uncertaiitty  in  Transition  Length 

Measuiements  of  tiansiticni  lengdi  at  a  given  heating  rate  were  conqtlicated  by 
the  fact  that  the  indication  of  transititm  start,  as  given  by  an  LC,  occurs  at  a  diffemit 
heating  rate  dian  does  dM  indicatitm  of  transition  end  (see  Sect.  3.1. 1.1).  Hieiefore, 
a  number  of  extrs^wlaticms  were  required  in  order  to  deteimine  the  locations  of 
transition  start  and  end  for  a  single  heating  rate.  The  procedure  is  described  in  detail 
in  Section  3. 1.1.2.  Briefly,  the  use  of  fluee  LC’s  provided  ReyntAds  numbers  at 
transition  end  for  three  different  heating  rates.  These  were  used  to  determine  the 
Reynolds  number  at  transition  end  for  zero  heating,  The  values  of  Rcr  and  Str 

for  each  liquid  crystal  were  dien  varied  individually  by  an  amount  equal  to  SRsr  and 
SStf  respectively  (as  shown  in  Table  B.2)  and  the  effect  on  was  observed.  This 
provided  a  measure  of  dRer^Rtr  and  dRer^dStj  such  that  die  uncertainty  in 
could  be  found  using  Eqns.  B.3  and  B.4.  Given  the  heating  rate  at  transition  start  and 
Rct^,  Etpi.  3.10  could  then  be  used  to  determine  the  location  of  sy  at  die  same  headng 
rate  as  x,.  The  diffetence  between  Xj  and  x,  is  the  transidon  lengdi,  T|,  at  a  given 
heating  rate.  Transition  loigth  was  normalised  with  die  distance  to  transidon  start. 

Uncertainty  in  the  correladon  whidi  quantifies  the  effect  of  headng  cm 
transidon  (Etpi  3.10)  was  estimated  based  on  a  95%  confidence  interval  constracted 
with  die  aid  of  Student’s  t-distribution.  Figure  3.14  shows  Ecpi.  3.10  superinqiosed 
on  all  experimental  measurements  of  transidon  Reynolds  number  versus  heating  rate 
at  zero  pressure  gradiem.  The  vertical  di^renoe  between  die  correladon  and  the 
experimental  data  was  calculated  for  each  poim  and  the  resulting  data  set  was  assumed 
to  be  normally  distriboted.  Student’s  t-distribudon  (with  41  degrees  of  freedom)  could 


then  be  used  to  constmct  a  95%  confidence  interval  for  the  correlation. 


Fig.  B.l  is  a  schematic  illustration  of  the  means  by  which  uncertainties 
conqwund  when  this  technique  is  applied.  Tte  numerical  uncertainties  included  in 
this  figure  represent  those  for  a  speed  of  22.3  m/s  with  a  heating  rate  of  St  =  9.23  x 
lO*^  (that  at  which  the  37.1°C  LC  indicated  the  location  of  x,).  The  relative 
uncertainties  in  the  calculated  values  of  normalised  transition  lengdis  for  the  23.7, 29.7 
and  37.1°C  LC’s  at  this  condition  were  ±  24.0,  22.0  and  21.9%  respectively.  Bars 
representing  these  uncertainties  are  superimposed  on  measurentents  of  normalised 
transition  length  in  Figure  3.17. 

Figure  3.14  illustrates  uncertainties  in  plots  of  Re/Re^  v.  5r  at  the  same 
conditions  as  discussed  above.  For  an  IXZ  mixtuie  with  response  tenqreratuies  of  23.7, 
29.7  and  37.1'’C,  unceitainties  in  Rej/RCr^  arc,  respectively,  6.4,  6.5  and  6.9%. 
Values  of  have  been  provided  in  Table  B.2. 

The  relative  uncertainty  in  u  is  slightly  larger  in  favourable  pressure  gradients 
since  local  velocity  was  not  measured  directly.  However,  as  fire  distance  to  transition 
start  also  increases  in  favourable  pressure  gradierrts  (transition  is  delayed),  die  smaller 
relative  uncertainty  in  x,  more  than  compensates  fer  the  additional  uncertainty  in  u. 
Therefore,  the  uncertainties  in  Re,  St,  RetRe^  and  T\lx,  discussed  here  are  a  worst-case 
illustration. 
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Appendix  C:  Estimates  of  Longitudinal  Conduction 

Estimates  of  longitudinal  conduction  due  to  suream-wise  suiface  tenqteratuie 
variations  on  the  heated  plate  ate  provided  in  this  qjpendix  for  both  the  aluminium 
and  carbon  heaters. 

Consider  an  mrinitesimal  element  of  a  heated  flat  plate  (representing  both 
heater  and  substrate)  of  unit  width  as  shown  in  Pig  C.l.  Assume  that  the  plate  is  in 
steady-state  and  that  die  longitudinal  ten^ratuie  distribution  is  uniform  in  die  span- 
wise  diiection.  Also  assume  that,  at  a  given  stream-wise  location,  the  temperature 
within  the  element  does  not  vary  with  y  (i.e.  T,  »  T/jc)).  The  heat  generated  within 
the  element  is 

=  f*  0  doc  ,  (C-1) 

where  (  is  the  current  and  9  is  the  local  lesistance/square  of  the  heater.  Heat  is 
conducted  into  the  left  face  and  out  the  right  face  of  the  element.  Heat  is  also 
convected  out  the  ui^r  suiface.  The  lower  suiface  is  assumed  to  be  perfectly 
insulated. 

d=d  ~d  j  1  (C,2) 

’f  ’C4M1V  'leone, 0*1  *tenJ,iH 

According  to  Fourier’s  law,  the  heat  flowing  into  the  left  face  of  the  lattice  can 
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be  expressed  as 


dr 


Similady,  conduction  out  the  right  side  of  the  lattice  is 


dr  d^r 

=  -fct  - ^dx 

.  d*  dx* 


Convection  may  be  represented  as 


4.,,,  =  hdxdT  =  hdxiT,-T^) 


Replacing  Eqn.  C.2  with  Eqns.  C.l,  C.3  and  C.4, 


,  dr  d^r  dr 

0  dx  =  ft  dx  (r  -  r)  -  ifc  f  — i  + - i  dx  +  ft  t  — i. 

*  *  dx  dx^  dx 


Sinq)lifying, 


i*e  =  ft(r  -  r.)  -kt — 

*  dx* 


Thus,  a  ratio  of  (2)  to  (1)  should  provide  a  measure  of  the  relative  in^rtance  of 
conduction  for  these  experiroeiits.  A  negative  value  of  (2)  represents  a  net  heat  loss 
from  the  element  under  consideration  due  to  conduction,  while  a  positive  value  of  (2) 
indicates  a  net  gain  in  heat.  Given  the  high  surface  temperature  characteristic  of 
transition  start,  negative  values  of  (2)  are  expected  in  diis  regitm.  Conversely,  (2) 
should  be  positive  near  transition  end. 


Conduction  will  be  largest  (in  relation  to  convection)  in  transitional  boundary 
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layen  at  low  speed.  At  reduced  flow  speeds,  convective  heat  transfer  is  relatively 
low,  and  temperature  gradients  are  highest  in  transitional  boundaiy  layers.  Figure  C.2 
shows  the  stream-wise  surface  temperature  distribution  along  a  flat  plate  with  a 
constant  heat  flux  surface  and  a  transitional  boundary  layer  at  24.2  m/s  as  predicted 
by  STANS.  For  this  {uediction,  experimental  evidence  was  used  to  obtain  the  iiqnit 
parameters  of  x,  =  0.16  m,  X]-  =  0.22  m  and  9,  =  1060  W  m'^. 

Attempts  to  analytically  di^rentiate  a  polynomial  fitted  to  the  line  in  Fig.  C.2 
did  not  produce  consistent  results.  Therefore,  a  numerical  second  derivative  of  die 
temperature  distribution,  (fjydx^,  was  determined  from  the  data  .  Results  are  shown 
in  Figure  C.3.  Conductive  losses  are  highest  slightly  downstream  of  transition  start, 
while  gains  due  to  conduction  are  largest  at  approximately  the  location  of  transition 
end. 


Using  Eqn.  C.7,  longitudinal  conductive  losses  were  found  to  be  a  maximum 
of  13.0%  of  convection  for  the  carbon  plate.  Maximum  conductive  gains  amounted 
to  12.6%  of  convectitHi.  The  measured  thickness  of  the  carbon  heater  (r)  was  3  mm, 
and  the  estimated  thermal  conductivity  of  the  carbon  fibres  was  6.7  Wm  '  -K  '  (see 
SecticHi  2.2,2.2).  Note  that  die  relative  importance  of  longitudinal  ctmduction  should 
be  much  less  in  the  nacelle  fli^t  tests  where  die  convective  heat  transfer  coefficient 
will  be  higher,  and  the  thermal  ccmductivity  of  the  carbon-fibre  matrix  reduced. 

Relative  longitudinal  conduction  in  the  perspex  plate  was  ccMisiderably  smaller 
than  in  the  carbon  plate  for  the  same  conditions.  Maximal  values  of  and 
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^ccmdjA  ■'*'***  0.8  and  0.7%  of  convective  heat  transfer  respectively.  The  effective 
thickness  of  the  perspex  plate  for  these  calculations  was  6  mm.  The  thermal 
conductivity,  k,  of  the  perspex  was  assumed  to  be  0.2  Wm'*'K'‘,  as  listed  by 
TOULOUKIAN.  et  al.  (1970). 

Conduction  does  not  affect  the  location  of  transition  as  indicated  by  the 
temperature-sensitive  LC’s  in  these  experiments.  Nevertheless,  conductive  Ic'ses  (or 
gains)  will  introduce  errors  into  local  heat  transfer  nteasurements.  Specifically, 
measurements  of  heat  transfer  coefficients  are  expected  to  be  spuriously  high  in  areas 
where  heat  is  being  lost  through  conduction  (transition  start).  The  reverse  is  true  in 
areas  where  heat  is  gained  via  conduction  (transition  end).  Since  the  adiabatic 
Reynolds  number  is  determined  with  the  aid  of  measurements  of  heat  transfer  at 
transition  start  and  end,  its  accurate  may  suffer  as  a  result  of  conduction.  A  lj% 
error  in  riKasutements  of  convective  heat  transfer  is  estimated  to  result  in  an  error  of 
6.2%  in  adiabatic  Reynolds  number  in  these  tests. 
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Figure  1.1  Extended  laminar  flow  due 
nacelle,  after  Hastings,  et  al.  (1987). 
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FIgvc  t2  Savings  due  to  1%  lednction  of  fuel  burned,  after  Szodruch  (1991). 


[ICTl 


Figure  Curve  of  neutral  stability  for  a  boundary  on  a  flat  plate  at  zero  incidence 
(Blasius  profile).  Figure  according  to  SchUchting  (1979);  theory  by  Tollmicn  (1929). 
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F%ve  lii  Turbulent  qxrt  rrutiated  by  electric  qraik  between  needle  electrode 
surface,  ato  Sdurbauer  and  Klebanoff  (1955). 
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Figure  1.7  Transition  via  the  T-S  path— &om  White  (1991) 


Figure  1.8  Three  dimensi<mal  development  of  a  single  wage,  depicted  in  consecutive 
wave  periods— from  Thomson  (1988) 


Figure  1.9  Three  dimensional  development  of  sueam-wise  vortices  downstream  of 
a  transverse  ti^  wire  in  a  laminar  boundary  layer.  Flow  visualisation  by  dye  in  water, 
after  Hama  et  al.  (1957). 
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Figure  1.19  Schematic  form  of  stream-wise  vortex  system,  after  Stuart  (1965). 


Figure  1.11  Instantaneous  velocity  profiles  across  boundary  layer  and  contour  of 
high-shear  stress  layer.  Measurements  of  Kovasznay  et  al.  (1962). 


Figure  1.12  Effect  of  surface  curvature  on  transition  Reynolds  number  in  zero 
pressure  gradient  flow,  after  Tani  (1%9).  Negative  values  of  S*/r  represent  convex 
curvature. 


Figure  1.13  Transition  Reynolds  number  as  a  functkm  of  fiee-stieam  tuibulenoe. 
Bjqwrimental  results  by  Diyden  (1959).  Line  represents  empirical  correlation  of  Mack 
(1977).  Figure  taken  fiom  Young  (1989). 


*  SchutMuer«ndSkxamsttd(I948) 

4  Hislop(1940) 

4  Abu  GhAniuiin  and  Shaw  (1980) 
a  Brown  and  Dunon  (1972) 

X  Manin.BrDwn  and  Gamtt  (1978) 

«  Dryden(l936) 

r  Costelow  and  Ramachandran  ( 1933) 

♦  Van  det  Hegje-Zijncn  (1924). 


Figure  1.14  Effect  of  hubulence  on  transition  Reynolds  number.  Correlation  from 
Govindarajan  and  Narasimha  (1991).  Note  expanded  scale  at  low  q. 


Heated  Surface 


Inflected  U(y) 
Transition  Hastened 


Favourable  Pressure  Gradient 


’Fuller’  U(y) 
Transition  Delayed 


Figure  1.15  Effects  of  beating  and  favourable  pressure  gradient  on  transition  in  an 
air  boundary  layer. 
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Figure  1.16  Experimentally  observed  influence  of  heating  on  transition  Reynolds 
numbers,  after  Lauchle  and  Gurney  (1984).  Axisymmetric  body  in  water. 
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Figure  1.17  Flat  plate  boundary  layer  velocity  profiles  in  the  transition  region,  as 
measured  by  Schubauer  and  Klebanoff  (1955). 
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Figure  1.18  Change  in  the  sh^  factor  H,j  of  the  boundary  layer  for  a  flat  plate  in 
the  transition  region  as  measured  by  Schubauer  and  Klebanoff  (1955)  and  taken  from 
SchUchting  (1979). 
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Figure  1.19  Transition  from  laminar  to  tuibulent  flow  as  indicated  by  successive  thin 
films,  after  Qaik,  et  al.  (1991).  Intermittency  «  0  at  flrst  gage  and  «>  1  at  last  gage. 


Figure  2«1  Low  tpetd  wind  tunnel  and  gcsieral  experimental  qiparatus. 


Figure  U1  rwamie  Fairing. 


a.  Typical  heat  transfer  distribution, 


b.  Low  heating  showing  transition  start. 


c.  High  heating  showing  transition  end. 


Figure  The  use  of  tempentuie  sensitive  li^iid  oystal  for  transition  detection. 


Figure  2.4  Flat  plate  and  aluminium  heater  used  for  heat  tiansfin'  expeiiments  (flow 
right  to  left,  600  mm  ruler  shown). 


Figure  2£  Electrical  contact  to  aluminium  heater. 


Figwc  2j6  Liquid  crysul/alominium  heatei  padcaging  anangement. 


F^«re  2S  Measured  vs.  actual  firee-stream  vdocity. 


Figure  2.9  Carbon  heater  mounted  in  perspex  plate  (iqiper  surface). 
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Figure  2.11  Carbon  heater  mounted  in  perspex  plate  Oower  surface).  Note; 
styrofoam  insulation  has  been  removed. 


FIgvc  2.12  Cross  sectioa  of  electrical  contact  on  carbon  heater. 


(a)  B^bre  soaUng  exposed  unidirectional  carbon  fibres  in  conducting  paint. 


(b)  Afier  soaHng  etq^osed  unidirectional  carbon  fibres  m  conducting  paint. 
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Figure  2.14  Vahatkm  of  diennal  conductivity  with  temperatuic  for  imiHimrrinnai 
caibon  fibres,  after  Rhodes  and  Moses  (1991). 


h(WmlC)  h(Wm^iC‘) 


n - 1 - 1  I  I 


KEY 

O  ?5.1  deg  cryslol 

*  33.8  deg  crystol 

•  37.4  deg 


'  ■  ■  ■  ' 


99  9.1  9.2  9.S  9.4  9.S  9.6  0.7  0.8  0.9 


KEr 

o  25.1  deg  crystal 
+  33.8  deg  crystol 
•  37.4  deg  crystal 


9.9  O.I  0.2  e.J  0.4  0.$  0.6  0.7  0  9  0.9 


Distance  along  plate  (m) 


Distance  along  plate  (m) 


(a)  242  misec 


(b)  28.9  misec 


1)9 

129 

♦  "1 

.  •  '■Ill - r 

. 

Ht 

. 

101 

- 

» 

4  •  • 

- 

81 

o  *  o* 

^9 

h  •'  o 

••  o 

60 

t  •, 

- 

SO 

40 

*•  V' 

KEY 

o  25  1  dCQ  crgstol 

♦  33.8  deq  crnistol 

•  37.4  deg  crystal 

—  »  » 

1  A  1  i  1  1  i 

9  9  0  1  0  2  O.S  9.4  O-S  9.6  9.7  0.8  0.9 


•  »  < 

.» 

o  o 

4  )(• 


KEY 

o  24.1  deg  crystal 
*  34.1  deg  crystal 
^  39.2  deg  crystal 


9.0  0.1  0.2  0  }  0.4  O.S  0.6  9.7  0.8  0.9 


Distance  along  plate  (m) 


Distance  along  plate  (m) 


(c)  32.2  m/sec 


(d)  34.8  misec 


Figure  3  J  Heat  transfn  distribution  on  heated  plate  (aluminiuin  heater). 


Figure  33  Influence  of  heating  on  transition  location.  □  24.2  m/s,  028.9  m/sec, 
A  32.2  m/sec,  034.8  m/sec.  Open  symbols— transition  start.  Oosed  symbols- 
transition  end.  Solid  line— least  squares  fit  to  data. 


Distance  along  plate  (m)  distance  along  plate 

(a)  /Nominal  speed  =  22  m/sec.  (f,)  Nominal  speed  =  28  mJsec. 
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(c)  Nominal  speed  =  35  m/sec. 


Figure  Favourable  pressure  gradient  on  flat  plate. 


Figure  3.6  Repeatability  of  favourable  pressure  gradients.  Nominal  speed  s  35 
m/sec.  Broken  and  solid  lines  represent  different  tests. 


Figure  3.7  Repeatability  of  favourable  piessuie  gradients.  Nominal  speed  =  28 
m/sec.  Broken  and  solid  lines  represent  diffeient  tests. 
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Variation  in  adiabatic  transition  location  with  pressure  gradient. 
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Figure  3.9  Influence  of  combined  heating  and  favourable  pressure  gradient  on 
transition  location.  Transition  start,  conditimi  1  (nominal  oncoming  speed  =  22 
m/sec). 


Figure  3.10  Influence  of  combined  heating  and  favouraUe  pressure  gradient  on 
transition  location.  Transition  start,  craidition  2  (nominal  oncoming  speed  =  28 
m^ec). 
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IFIgure  3JL1  Influence  of  combined  beating  and  favourable  pressure  gradient  on 
transition  locatitni.  Transition  end,  condition  2  (nominal  oncoming  speed  =  28  m/sec). 
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F%vre  3.12  Influence  of  combined  heating  and  favourable  pressure  gradient  on 
transition  location.  Tianntion  stait,  cooditwn  3  (nominal  osKoming  speed  »  35 
m/sec). 
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Fignrc  3.15  Influence  of  heating  and  piessuie  gradient  on  transition  location. 
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Figure  3.17  Influence  of 
length. 
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Figure  3.19  Heat  transfo'  distribution  on  caibon  heatra.  Flow  speed  =  34.3  m/sec. 


Figure  3.21  Comparison  of  heat  transfer  distributions  for  trqtped  and  untripped 
boundary  layers  on  carbon  plate.  Open  symbols— untrqtped.  Oosed  symbols— tripped. 
Flow  speed  =  28.9  m/sec. 


Figure  3.22  Comparison  of  heat  transfer  distributions  for  tripped  and  untripoed 
boundary  layers  on  carbon  plate.  Open  symbols— untripped.  Qosed  symbols— tr^ped. 
Flow  =  7.9  m/sec. 


Figure  3J3  Contours  of  constant  temperattu'e  on  carbon-fibre  heater.  Turbulent 
boundary  layer. 


Figure  3.24  Transitional  boundary  layer  on  carbon-fibre  heater  as  indicated  by 
temperature  sensitive  liquid  crystals. 
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Figure  3.25  Heat  transfer  distribution  on  carbon  plate.  Transition  boundary  layer. 
Transition  stan  for  prediction  upstream  of  the  heated  section.  Flow  q>eed  =  28.9  m/s. 


Figure  3M  Turbulent  wedge  cm  heated  carbon  (date. 
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Fisure  3^7  Time  required  for  carbon  plate  to  teach  thermal  stability.  Flow  speed 
=  28.9  m/sec. 
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Figure  3.28  Time  required  for  carbon  plate  to  teach  thermal  stability.  Flow  ^ed 
=  7.9  m/sec. 


Figure  3J9  Curved  caibon-fibie  heater. 


Figure  3  J1  Transhioaal  boundary  layer  on  curved  carbon-fibre  heater. 
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Figure  3J3  Results  from  two  tuns  at  36^  m/sec  and  STANS  prediction  (with 
specified  transition  location)  used  to  calibrate  the  liquid  crystal.  This  calibration  is 
then  used  to  determine  the  shear  stress  for  two  runs  at  33.S  m/sec. 
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Repeatability  of  shear  sensitive  liquid  crystal  tests.  Nominal  q;>eed  = 


5.0 

0.5 

- 

— 1 - i  1  1  \ 

1 

0.0 

- 

- 

3.5 

- 

- 

3.0 

o 

•  O  O 

- 

3.5 

X  o  .  .  _ 

. 

3.0 

- 

Op  X  o  • 

?  -  o 

- 

°  Tq  =  20.0  deq.  C 

I.O 

*  Tq  19.3  deq.  C 
+  Tq  =  18.8  deq.  C 

0.5 

- 

o  ig  =  18.0  deq.  C 
^  Tg  =  16.0  deg.  C 

- 

D.O 

1  1  1  1  1  i 

1 

C.C  C.l  0.3  0.3  o,*.  O.S  0.6  0.7  0.8  C.6 


Distance  along  plate  (m) 


Figore  3J5  Repeatability  of  shear  sensitive  liquid  ciystal  tests.  Nominal  qoeed 
33.6  m/sec. 


Flgarc  3M  ‘Hirtxileiit  wedge  in  laminar  boundary  layer  detected  with  drear  sensitive 
liquid  crystal.  Flow  tpced  *  36  j  m/sec. 
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Figure  339  T^iirbulent  wedge  produced  by  three-dimensional  surface  excrescence, 
after  Schubauer  and  Klebanoff  (1953). 
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Figure  3.40  Contours  of  heat  transfer  in  a  turbulent  wedge,  after  Clarir  and  Jones 
(1991). 


Figure  3.41  Stream-wise  vortices  downstream  of  a  spherical  protuberance  on  a  flat 
plate,  after  Hama  et  al.  (1957).  Flow  visualisation  by  sublimating  acoraphthene. 
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Figure  A.l  Transitional  boundary  layer  on  a  uniformly  heated  flat  plate  in  an 
incompressible  fluid. 


